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Two-sided bounds for degenerate processes with densities 
supported in subsets of M. N 

Chiara Cinti? Stephane Menozzi "^and Sergio Polidoro* 



O. 

fN| , Abstract: We obtain two-sided bounds for the density of stochastic processes satisfying a weak 

q ' Hormander condition. In particular we consider the cases when the support of the density is 

not the whole space and when the density has various asymptotic regimes depending on the 

f— 1 [ starting/final points considered (which are as well related to the number of brackets needed 

to span the space in Hormander 's theorem). The proofs of our lower bounds are based on 

Harnack inequalities for positive solutions of PDEs whereas the upper bounds are derived 

from the probabilistic representation of the density given by the Malliavin calculus. 
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1 Introduction 

We present a methodology to derive two-sided bounds for the density of some Revalued 
degenerate processes of the form 

ov _ 

en ' 

O 

CN ' where the (^)ie[o,n] are smooth vector fields defined on R^, {(W^)t>o)ieU,nl stand for 

standard monodimensional independent Brownian motions defined on a filtered probability 
space (fl, J£ ', (^t)t>o, P) satisfying the usual conditions. Also o dWt denotes the Stratonovitch 
integral. The above stochastic differential equation is associated to the Kolmogorov operator 

a" 



l [ Yi(X s ) o dWt + / Y (X s )ds (1.1) 
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^=\Y. Y ? + Z > Z = Y -d t . (1.2) 
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We assume that the Hormander condition holds: 
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[H] Rank(Lie{Yi,--- ,Y n ,Z}(x)) =N + 1, Vx G 
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We will particularly focus on processes satisfying a weak Hormander condition, that is 
Rank(Lie{Y"i, • • • ,Y n , —dt}(x)) < N + 1, Vx G M. N . This means that the first order vector 
field Yq (or equivalently the drift term of the SDE) is needed to span all the directions. 

As leading examples we have in mind processes of the form 

X\ = Xl + Wl, Vi G [l,n], A7+ 1 = x n+1 + f \Xl' n \ k ds, (1.3) 

Jo 

where X s ' n = (Xj, • • • , X™) (and correspondingly for every x G R n+1 , x\^ n := (x\, • • • , x n )), 
k is any even positive integer and |.| denotes the Euclidean norm of W 1 . Note that we only 
consider even exponents in f|l .3j) in order to keep Yq smooth. Our approach also applies to 

,-t n 
X$= Xi + Wi, ViG[l,n], X? +1 = x n+1 + / y(Xl) k ds, (1.4) 
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for any given positive integer k. 

It is easily seen that the above class of processes satisfies the weak Hormander condition. 
Also for equation (|1.3|) . the density p(t, x, .) of X t is supported on W 1 x (x n+ i, +oo) for any 
t > 0. Analogously, for equation (|1.4p . the support of p(t,x,.) is R n+1 when k is odd and 
M. n x (x n +i, +oo) when k is even. 

Let us now briefly recall some known results concerning these two examples. First of all, 
for k = 1, equation (|1.4p defines a Gaussian process. The explicit expression of the density 
goes back to Kolmogorov [25J and writes for all t > 0, x,£ G M n+1 : 

V3 / I 1 |£i )Tl — Xl, n | „|?n+l — %n+l 2 *" 



pit(<,x,o = (27r) ^^ ex p^- ji ' 1 ' + 3 p — = ^ 1 -a- 5 ) 

We already observe the two time scales associated respectively to the Brownian motion (of 
order t 1 ' 2 ) and to its integral (of order t 3 ' 2 ) which give the global diagonal decay of order 
£n/2+3/2_ rp^g a( j^iti ona i term Xl 2 t in the above estimate is due to the transport of the 
initial condition by the unbounded drift. We also refer to the works of Cinti and Polidoro 
|17| and Delarue and Menozzi |19j for similar estimates in the more general framework of 
variable coefficients, including non linear drift terms with linear growth. 

For equation (|1.3|) and k = 2, n = 1, a representation of the density of Xf has been 
obtained from the seminal works of Kac on the Laplace transform of the integral of the 
square of the Brownian motion [23 . We can refer to the monograph of Borodin and Salminen 
[10J for an explicit expression in terms of special functions. We can also mention the work of 
Tolmatz [39] concerning the distribution function of the square of the Brownian bridge already 
characterized in the early work of Smirnov [36]. Anyhow, all these explicit representations 
are very much linked to Liouville type problems and this approach can hardly be extended 
to higher dimensions for the underlying Brownian motion. Also, it seems difficult from the 
expressions of JTO] to derive explicit quantitative bounds on the density. 



Some related examples have been addressed by Ben Arous and Leandre |5j who obtained 
asymptotic expansions for the density on the diagonal for the process X\ = x\ + W t , X 2 = 
%2 + fo{Xl) m dWs + f (X],) k ds. Various asymptotic regimes are deduced depending on m 
and k. Anyhow, the strong Hdrmander condition is really required in their approach, i.e. the 
stochastic integral is needed in X 2 . 

From the applicative point of view, equations with quadratic growth naturally appear in 
some turbulence models, see e.g. the chapter concerning the dyadic model in Flandoli |20\ 
This model is derived from the formulation of the Euler equations on the torus in Fourier 
series after a simplification consisting in considering a nearest neighbour interaction in the 
wave space. This operation leads to consider an infinite system of differential equations 
whose coefficients have quadratic growth. In order to obtain some uniqueness properties, a 
Brownian noise is usually added on each component. In the current work, we investigate from 
a quantitative viewpoint what can be said for a drastic reduction of this simplified model, 
that is when considering 2 equations only, when the noise only acts on one component and 
is transmitted through the system thanks to the (weak) Hdrmander condition. 

Our approach to derive two-sided estimates for the above examples is the following. The 
lower bounds are obtained using local Harnack estimates for positive solutions of J£u = 
with Jz? defined in (jl.2p . Once the Harnack inequality is established, the lower bound for 
pit, x, £) is derived applying it recursively along a suitable path joining x to £ in time t. The 
set of points of the path to which the Harnack inequality is applied is commonly called a 
Harnack chain. For k = 1 in (|1.4|) the path can be chosen as the solution to the deterministic 
controllability problem associated to (|1.4p . that is taking the points of the Harnack chain 
along the path 7 where 

n 

7-(s) = uji(s), Mi e [l,n], 7n+i( s ) = J^ 7 *^' 7 ^ = x ' 7< ^ = ^ 

i=i 

and uj : L 2 ([0, £]) — > W 1 achieves the minimum of L \u(s)\ 2 ds, see e.g. Boscain and Polidoro 
|llj . Carciola et al. [33] and Delarue and Menozzi |19j . 

In the more general case k > 1 it is known that uniqueness fails for the associated control 
problem, i.e. when 7' +1 (s) = X^=i(7i( s )) m the above equation (see e.g. Trelat jlTFj). 
Therefore, there is not a single natural choice for the path 7. Actually, we will consider 
suitable paths in order to derive homogeneous two-sided bounds. After the statement of our 
main results, we will see in Remark 12.31 that the paths we consider allow to obtain a cost 
similar to the one found in [3D] for the abnormal extremals of the value function associated 
to the control problem. 

Anyhow, the crucial point in this approach is to obtain a Harnack inequality invariant 
w.r.t. scale and translation. Introducing for all (m,x) G N* x W l+l the space V m (x) := 

{(0 / il)tie[l,n]> (Kn^ 2 ]( a; ))(Ji,i2)G[0,n] 2 r • • > (\Yii, [ Y 'i 2 r-- , [ Y i m -n Y i m }]}( x ))(h,-,im)&lO,n} m )}> 

the above invariance properties imply that dim(Span{V^,(x)}) does not depend on x for any 
m. This property fails for k > 1 since we need exactly k brackets to span the space at 
x = (Oi jn ,x n +i) and exactly one bracket elsewhere. Hence, we need to consider a lifting 
procedure of Jzf in (|1.2|) introduced by Rotschild and Stein |B_5j (see also Bonfiglioli and 
Lanconelli |6j). Our strategy then consists in obtaining an invariant Harnack inequality 



for the lifted operator 5£ . We then conclude applying the previous Harnack inequality to 
££ -harmonic functions (which are also «Sf-harmonic) . A first attempt to achieve the whole 
procedure to derive a lower bound for (|l,4p and odd k can be found in Cinti and Polidoro 

[TBI. 



Concerning the upper bounds, we rely on the representation of the density of p obtained by 
the Malliavin calculus. We refer to Nualart |33| for a comprehensive treatment of this subject. 
The main issues then consist in controlling the tails of the random variables at hand and the 
IP norm of the Malliavin covariance matrix for p > 1. The tails can be controlled thanks to 
some fine properties of the Brownian motion or bridge and its local time. The behavior of the 
Malliavin covariance matrix has to be carefully analyzed introducing a dichotomy between 
the case for which the final and starting points of the Brownian motion in (|1.3p - (|1.4p are close 
to zero w.r.t. the characteristic time-scale, i.e. |a?i n | V |£i n| 5; Kt 1 ' 2 for a given K > 0, which 
means that the non-degenerate component is in diagonal regime, and the complementary set. 
In the first case, we will see that the characteristic time scales of the system (|1.3p . (|l-4p and 
the probabilistic approach to the proof of Hormander theorem, see e.g. Norris [3JJ will lead 
to the expected bound on the Malliavin covariance matrix whereas in the second case a more 
subtle analysis is required in order to derive a diagonal behavior of the density similar to 
the Gaussian case (|1.5p . Intuitively, when the magnitude of either the starting or the final 
point of the Brownian motion is above the characteristic time-scale, then only one bracket is 
needed to span the space and the Gaussian regime prevails in small time. 

Note that our procedure can be split in two steps. In the first one, purely PDEs methods 
provide us with lower bounds of the density p. In this part useful information about its 
asymptotic behavior in various regimes are obtained by elementary arguments. Once the 
lower bounds have been established, we rely on some ad hoc tools of the Malliavin calculus 
to prove the analogous upper bounds. However, aiming at improving the readability of our 
work, we reverse our exposition: we first prove the upper bounds, as well as the diagonal 
ones, by using probabilistic methods, then we prove the lower bounds by PDEs arguments. 

The article is organized as follows. We state our main results in Section [2j We then recall 
some basic facts of Malliavin calculus in Section and obtain the upper bounds as well as a 
diagonal lower bound in Gaussian regime in Section Hlln Section [5j we recall some aspects of 
abstract potential theory needed to derive the invariant Harnack inequality. We also give a 
geometric characterization of the set where the inequality holds. Section [6] is devoted to the 
proof of the lower bounds. 

2 Main Results 

Before giving the precise statement of our bounds for the the density p of X in (|1.3|) or (|1.4p . 
we give some remarks. In the sequel p(t, x, .) stands for the density of any stochastic process 
X at time t starting from x. It is well known that, if the vector fields Y\, . . . ,Y n (note that 
the drift term Yq does not appear) satisfy the Hormander condition, then the following two 
sided bound holds: 



Here and in the sequel, for measurable functions g : R + * x R n — >• R, h : R + * x R 2n , the above 
notation p(i, x,£) x -75— rexp(— /t(i, x,£)) means that there exists a constant C > 1 s.t. 

(7—1 (7 

■ exp(-C7i(t, x, £)) < P(*, a;, < —, r exp(-C _1 /i(t, x, £))• (2.2) 



S'(*,») ' " ' g(t,x) 



Moreover, in (|2.ip . dy denotes the Carnot-Caratheodory distance associated to Y\,...,Y n , 
and By(x,r) is the relevant metric ball, with center at x and radius r. On the other hand 
(|1.5p shows that, when the drift term lo is needed to check the Hormander condition, the 
density p of the process X doesn't satisfy ()2.ip . In this article we prove that, when considering 
processes (|1.3|) and (|1.4p with k > 1, different asymptotic behavior as \x\ — > +00 appear. 

To be more specific, we first remark that a behavior similar to (|1.5p can also be observed 
for equations (|1.3|) and (| 1 . 4p .Conditioning w.r.t. to the non degenerate component we get 

p(t,x,£) =p x i,n(t,xi, n ,£,i, n )p X n+i(t,x n+1 ,£, n+1 \X ' n = x ltn ,X t ' n = £ 1>n ), 



where Px x >™ (i> £i,raj £i,n) = (2 An/2 ex P( 2^ ) * s ^ e usua l Gaussian density, 

-l.n -v-X.n 



(2irt)' 1 
PXn+l(t,X n+ l,£ n+ l\XQ' n =Xi )Tl ,Xl' n = £l,n) =PY t {in+l — X n+ l), 



Y \Sl \^ Xl ' n + f&.« + w ^\ kds for C3D' 

1 " ' Uo EILi(¥^ + §& + W^) k ds for dHD, 

and (W s ' )ue[o,t] stands for the standard d-dimensional Brownian bridge on [0, t], i.e. starting 
and ending at 0. 

For the sake of simplicity, we next focus on the case n = 1 and k = 2 so that (|1.3|) and 
(|1.4p coincide. Moreover we assume x\ = £1. This leads to estimate the density of: 

Y t := I (xi + W^) 2 ds = tx\ + 2xi / W°'*ds + f {W^fds. (2.3) 

Jo Jo Jo 

Thus, when |xi| is sufficiently big w.r.t. the characteristic time scale t 1 ' 2 , the Gaussian 
random variable 

G := 2 Xl I'' W^ds ^ M (0,^\ Xl \ 2 

dominates in terms of fluctuation order w.r.t. the other random contribution whose variance 
behaves as 0(t A ) q 

If we additionally assume that |^2 — X2 — tx 2 \ < C\x\\t 3 ' 2 , for some constant C := C(n = 
1, k = 2) to be specified later on, that is the deviation from the deterministic system deriving 



lr The previous identity in law is derived from Ito's formula and the differential dynamics of the Brownian 
bridge. Namely, /„' W°'*ds = {-(«- s)W°% + f*(t - s) (-^*> + dW s ) } <=* ft W°'*da = \ f*(t - 
s)dW s . 



from (|1.3p . obtained dropping the Brownian contribution, has the same order as the standard 
deviation of G, we actually find: 



1+3— fi , i UJ2+3 

X\\t 2 VI'' l X U L 

When |^2 — X2~ tx\\ > C|xi|i 3 ' 2 , that is when the deviation from the deterministic system 
exceeds a certain constant times the standard deviation, the term J (W S ' ) 2 ds in (|2.3|) is not 
negligeable any more and we obtain the following heavy-tailed estimate: 

^.^^{M^ + Mi^M 

The diagonal contribution of the degenerate component corresponds to the intrinsic scale of 
order t 2 of the term L (W s ' ) 2 ds. In particular, if x\, n = 0i jn this is the only random variable 
involved. The off-diagonal bound can be explained by the fact that f Q (W s ' ) 2 ds belongs to 
the Wiener chaos of order 2. The tails of the distribution function for such random variables 
can be characterized, see e.g. Janson |22j, and are homogeneous to the non Gaussian term 
in the above estimate. 

Observe also that the density p is supported on the half space {£ £ M 2 : £2 > x 2 }- We 
obtain as well an asymptotic behavior for the density close to the boundary. Precisely, for 
< £2 — X2 sufficiently small w.r.t. to the characteristic time-scale t 2 of L (W s ' ) 2 ds, that is 
when the deviations of the degenerate component have the same magnitude as those of the 
highest order random contribution, then 

( + t\~ X ( J >i| 4 + I6| 4 , t 2 



1 1/2+2 ^\, I 6-^2 (6-X 2 ) 

We summarize the above remarks with the assertion that processes of the form (|1.3|) or 
(|1.4|) do not have a single regime for k > l.The precise statements of the previous density 
bounds are formulated for general n and k in the following Theorem 12.11 



Theorem 2.1 Let x = (xi >n ,x n+ i) G R n+1 , and £ = (£i, n >£n+l) G K n x (x„ +1 , +00) for k 
even, and £ G M. n+1 for k odd, be given. Define 



\xi, n \ k + |£l,n| fe , for ([L3D, 



i) Assume Nlv^^rte^TT ~ > C w/iere c := c(fc) = 2 + ^ and C is fixed. 
T/ien i/iere exists a constant C\ := C\{n, k, C) > 1 s.i. /or every t > 0, 

^^exp(-C 1 /(t,x,e,^))<p(t,x,0<^^exp(-Cr 1 /(t,x,e,|-^)),(2.4) 
Vc€R+ !(*,*, &c) := ^ 7 1 '"' 2 + lg " +1 ~ Xn+1 ~ l f 2 ^ il ' n)t? ' k • 



ii) Assume ^s^u* 1 pfe^jfey^ < C (with c,C as in point %)) and |xi in |V|6,n|A 1/2 > 

K, with K sufficiently large. Then, there exists C2 '■= Ci (n, k, K, C) > 1 s.t. for every 
t>0: 

C 2 l e W (-C2l(t,x,Q) n(frf)< g 2 exp(-C 2 " 1 J(t^,Q) 

n + 3 — r\ ' "^'v — . 1 , _, ra + 3 ? \ / 



(I^Lnl*- 1 + l^nl*- 1 )*"* 2 " (lxi.nl*" 1 + l^nl^ 1 )^ 

... ^ _ l^l,n ~ Xl,n| 2 |gn+l ~ X n+ l - ^(xi, n , £l,n)*| 2 

l ' X ' U ' t + (|xi, B |(*-i) + |ft, B |(»-i))V ' 

mj For t > 0, assume |£ n +i — x n +i| < Kt 1+k ' 2 for sufficiently small K. Then, there exists 
C3 := C^{n,k,K) > 1 s.f. we /iaue: 

C L eX p(_c 3 /(t,x,0) < p(t,x,0<^-exp(-C 3 - 1 /(t,x,0), 



£— + 1 £~ + 1 

| T , |2+fc i It, |2+fc fl+2/fc 

/(*,*,« := |X1 '"J + Igl>n + ^ (5^. (2.6) 

|4„+i - avj.il |Wi - x n+ i|^/' c 

Remark 2.2 /£ zs known that the densities of the stochastic systems (jl.3p and ()1.4[) are 
Fundamental Solutions of the Kolmogorov operators 

J? = ±A Xhn + \ Xl , n \ k d Xn+1 -d u (2.7) 

and 

n 

J? = -A xl>n +Y / xp Xn+1 -d u (2.8) 

i=i 
respectively. Then, from the PDEs point of view, Theorem \2.1\ provides us with estimates 
analogous to those due to Nash, Aronson and Serrin for uniformly parabolic operators. 

We next give some comments about our main result. As already pointed out, processes of 
the form (|1.3p or (jl.4|) do not have a single regime anymore for k > 1. Let us anyhow specify 
that when C~ l ^/i < \x{\ < Cy/i, \/i £ [l,n], C > 1, then expanding Yj as in (|2.3p . we find 
that all the terms have the same order and thus a global estimate of type (|2.4|) (resp. of type 
(|2.5|) ) holds for the upper bound (resp. lower bound) in both cases (|1.3p and (|1.4p . Observe 
also that in this case (|2.4|) and (|2.5|) give the same global diagonal decay of order £(«+ n )/ 2 +i. 

Remark 2.3 As already mentioned in the introduction, for k = 2,n = 1, we observe from 
(|2.6p that the off- diagonal bound is homogeneous to the asymptotic expansion of the value 
function associated to the control problem at its abnormal extremals, see Example 4-2 in \l~0j - 

The optimal cost is asymptotically equivalent to 44 1 when x = (0,0) as £ is close to (0,0). 

Remark 2.4 Fix |£ n+1 - x n+1 \ small, t G [K -1 |£„ + i - x n+1 \ 2 ~ £ , K\£ n+t - x n+ i\ 2 ~ £ ] for 
given K > l,e > 0. We then get from (|2.6p that there exist c := c(n,k),C := C(n,k,T) 
s.t. p(t,x,£) < Cexp(— c/|£ n +i — x n+ i| e ). This estimate can be compared to the exponential 
decay on the diagonal proved by Ben Arous and Leandre in Theorem 1.1]. 



3 A Glimpse of Malliavin Calculus 

3.1 Introduction 

Introduced at the end of the 70s by Malliavin, [30], [29], the stochastic calculus of variations, 
now known as Malliavin calculus, turned out to be a very fruitful tool. It allows to give 
probabilistic proofs of the celebrated Hormander theorem, see e.g. Stroock [37J or Norris 
[3T] . It also provides a quite natural way to derive density estimates for degenerate diffusion 
processes. The most striking achievement in this direction is the series of papers by Kusuoka 
and Stroock, |26j, [27], [28]. Anyhow, in those works the authors always considered "strong" 
Hormander conditions, that is the underlying space is assumed to be spanned by brackets 
involving only the vector fields of the diffusive part. For the examples (jl.3p . (|1.4p we consider, 
this condition is not fulfilled. Anyhow a careful analysis of the Malliavin covariance matrix 
will naturally lead to the upper bounds of Theorem 12.11 and also to a Gaussian lower bound, 
when the initial or final point of the non-degenerate component is "far" from zero w.r.t. 
the characteristic time scale on the compact sets of the underlying metric, see point ii) of 
Theorem 12.11 

We also point out that because of the non uniqueness associated to the deterministic 
control problem, the strategy of [19] relying on a stochastic control representation of the 
density breaks down. For the systems handled in [TO], we refer to Bally and Kohatsu-Higa for a 
Malliavin calculus approach [2 J . The Malliavin calculus remains the most robust probabilistic 
approach to density estimate in the degenerate setting. 

We now briefly state some facts and notations concerning the Malliavin calculus that 
are needed to prove our results. We refer to the monograph of Nualart [32], from which we 
borrow the notations, or Chapter 5 in Ikeda and Watanabe [21] . for further details. 

3.2 Operators of the Malliavin Calculus 

Let us consider an n-dimensional Brownian motion W on the filtered probability space 
(0,^",(^t) t >o,P) and a given T > 0. Define for h G L 2 (R+,R n ), W(h) = f^(h(s),dW s ). 
We denote by S the space of simple functionals of the Brownian motion W, that is the 
subspace of L 2 (Q,,^,¥) consisting of real valued random variables F having the form 

F = f(W(hi),--- ,W(hm)), 

for some m G N, hi G L 2 (R + , R n ), and where / : R m — > R stands for a smooth function with 
polynomial growth. 

Malliavin Derivative. 

For F G S, we define the Malliavin derivative (DtF) te iQ^] as the R n -dimensional (non 
adapted) process 



D t F = J29xj(W(h 1 ),--- ,W(h m ))hi(t). 



i=i 



For any q > 1, the operator D : S — > L q (Q,,L 2 (0,T)) is closable. We denote its domain by 
ED 1 ' 9 which is actually the completion of S w.r.t. the norm 



\F\\ l!q :={E[\F\i]+E[\DF\l 2{0jT) }] 



1/9 



Writing F)\F for the j component of DtF, we define the k order derivative as the random 
vector on [0, T] k x Q, with coordinates: 

D\yfF:=D\ k ---D\\F. 
ti, ■■• ,tk tfc ti 

We then denote by D ' 9 the completion of 5 w.r.t. the norm 

f N \ V« 

||F||jv,:= ^E[|FH + VE[|£> fe F| 9 , J I . 

II IMV.q | LI I j ^ LI 'LZ^o^fcyj 

Also, ED°° := n g >irij>ilIP'9. In the sequel we agree to denote for all q > 1, \\F\\ q := E[\F\ q ] 1 / q . 

Skorohod Integral. 

We denote by V the space of simple processes, that is the subspace of L 2 ([0, T]x!l,#x 
B([0,T]),dt <g> dP) consisting of W 1 valued processes processes (ut)te[o,T] that can be written 



u t 



^FiiWfa),--- ,W(h m ))hi(t), 



i=\ 



for some m G N, where the (-^Oieli.ml are smooth real valued functions with polynomial 
growth, V* G [l,m], /i* G L 2 ([0,T],IR™) so that in particular Fi(W(hi), ■ ■ ■ ,W(h m )) G S. 

Observe also that with previous definition of the Malliavin derivative for F G S we have 
(DsF) S £\Q t T] ^ T • F° r u £ V we define the Skorohod integral 

m , m 

8(u) :=^<U(W(^),--- ,W(h m ))W(hi) -Y.djFiiWihi),- ■ ■ ,W(h m )){h t ,h 3 ) L 2 {[m 

i=l *> j=l 

so that in particular 5(u) G S. The Skorohod integral is also closable. Its domain writes 
Dom(<5) := {u G L 2 ([0,T] x U) : 3(u n ) neP , u n L2([ ^i xn) Uj <5( Un ) L ^ F ■= S ( u )}. 

n n 

Ornstein Uhlenbeck operator. 

To state the main tool used in our proofs, i.e. the integration by parts formula in its whole 
generality, we need to introduce a last operator. Namely, the Ornstein-Uhlenbeck operator 
L which for F G S writes: 

LF := 5(DF) = (Vf(W(h)),W(h))-Tr(D 2 f(W(h))(h,h*) L 2 {0 ,T)), 
W{h) = (Wih),--- ,W(h m )). 



This operator is also closable and B°° is included in its domain Dom(L). 
Integration by parts. 



Proposition 3.1 (Integration by parts: first version) Let F G B 1,2 , u G Dom(5), then 
the following indentity holds: 

E[(DF,u) L 2 {[0>ri) ]=E[FS(u)], 

that is the Skorohod integral 5 is the adjoint of the Malliavin derivative D. As a consequence, 
for F,G G Dom(L) we have 

E[FLG] = E[F5(DG)} = E[(DF,DG) L 2 {[0jT]) ] = E[LFG], 

i.e. L is self-adjoint. 

These relations can be easily checked for F,G G 5, u G V, and extended to the indicated 
domains thanks to the closability. 

3.3 Chaos Decomposition 

Im(fm)--=m\ / ••• / f m (tl,'~ ,t m ) ® dW tm ® ■ ■ ■ ® dW tl . 

Jo Jo Jo 

In the above equation <g> denotes the tensor product and (dWt m ® • • • <8> dWt x ) G ((M n )® m )*. 
We now state a theorem that provides a decomposition of real-valued square-integrable 
random variables in terms of series of multiple integrals. 

Lemma 3.1 Let F be a real-valued random variable in L 2 (Q, J^,P). There exists a sequence 
(/mjmeN s.t. 

F=Y J Im(fm), (3.1) 

meN 

where for all mGN, f m is a symmetric function in L 2 ([0,T] m , (R n )® m ) and 



E[F \ — 2_^ m --||/m|li2([o,T] m ,(M n )® m ) ^ +00- 



m>0 



We refer to Theorem 1.1.2 in Nualart [32 J for a proof. 

Remark 3.2 We use the term chaos decomposition for the previous expansion because the 
multiple integral I m maps L 2 ([0,T] m , (R n )® m ) onto the Wiener chaos U m ■= {H m (W(h)), h G 
L 2 ([0,T],M n ), ||/i||l2qo,t],R") = 1}> where H m stands for the Hermite polynomial of degree m 
(see again Theorem 1.1.2 in [32]). The orthogonality of the Hermite polynomials yields the 
orthogonality of the Wiener chaos, i.e. E[XY] = 0, for (X, Y) G (T-L n ,T-L m ), n^m. 
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The computation of Malliavin derivatives is quite simple for multiple integrals. Indeed, 

Dt(I m (fm)) = ml ro -i(/ ra (i,.))el n . 
As a consequence, for a random variable F having a decomposition as in (|3.ip . we have 
that it belongs to B 1 ' 2 if and only if > J rnrn] -\\fm\\\2(\Q T ]- m mn)®™,) < +00 in which case 



m>l 
rT 



D t F = Y,m>i mI m-i{fm{t, • )) and E[/ \D t F\ 2 dt] = ^mm!||/ m || 2 2([ 0]T ]m( R n)®m). Iterating 

m>l 

the procedure, one gets F G O^' 2 «=>- Ylm=N (m%)\ ll/™lli,3([o,T]"s(R")®™) < +°° and 

+00 
D tu ... ttN F = Y^m(m-l)---(m-N + l)I m - N (f m (ti, ■ ■ ■ , t n , .)) G (R n ) 9N . 

m=N 

Therefore, when a random variable is smooth in the Malliavin sense, i.e. D°°, the Stroock 
formula, see [38], provides a representation for the functions (f m )meN i n the chaotic expansion 
in terms of Malliavin derivatives. 

Proposition 3.2 (Stroock's formula) Let F G B°°, then the explicit expression of the 
functions (/ m ) m >i in the chaotic expansion f)3.1|) of F writes: 

Vm G N, f m (tx, ■■■,t m )= E[D^... >tm F] G (R n f m . 

For square integrable process, a result analogous to Lemma 13. II also holds. 

Lemma 3.3 Let (ut)te[o,T) be an W l -valued process in L 2 ([0,T] x VL,& x B([Q,T]),dt <g> cflP). 
There exists a sequence of deterministic functions (g m )m£N* s.t. 

u t = ^2l m {gm+i(t, .)), (3.2) 

m>0 

where the square integrable kernels g m +i are defined on [0,T] m+1 with values in (]R n )®( m + 1 ) ) 
are symmetric in the last m variables and s.t. E m >o m! Hfi , m+ill!2([ ,T]™+i,(R™)®(™+i)) < + °°- 

We refer to Lemma 1.3.1 in |32] for a proof when n = 1. 

Also, the Skorohod integral of u G Dom(5) is quite direct to compute from its chaotic 
decomposition ()3.2p . Namely, 

<H n ) : = 5Z -Wi(Sm), 



m>0 



where g m (t,tx, ■ ■ ■ ,t m ) :— ^+i[9m(ti, • • • ,t m ,t) + Y^ILi 9m(ti, ■ ■ ■ ,U-i,t,ti+i, ■ ■ ■ ,t m ,ti 
the symmetrization of g m in [0, T] 



is 



m+1 

3.4 Representation of densities through Malliavin calculus 

For F = (Fi, ■ ■ ■ ,-F/v) G (D°°) , we define the Malliavin covariance matrix jp by 

lF j := (DF i ,DF^) L 2 i0>T) ,y(i,j) G [l,iVj 2 . 
Let us now introduce the non-degeneracy condition 
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[ND] We say that the random vector F = {F\, • • • , Fn) satisfies the non degeneracy condition 
if jf is a.s. invertible and det(7^) _1 G n q >iL q (Q). In the sequel, we denote the inverse of 
the Malliavin matrix by 

This non degeneracy condition guarantees the existence of a smooth density, i.e. C°° , for the 
random variable F, see e.g. Corollary 2.1.2 in [32] or Theorem 9.3 in [21] . 

The following Proposition will be crucial in the derivation of an explicit representation of 
the density. 

Proposition 3.3 (Second integration by parts) Let F = (i*i, • • • ,-F/v) G (B 00 )^ satisfy 
the nondegeneracy condition [ND] . Then, for all smooth function ip with polynomial growth, 
G G B°° and all multi-index a, 

K[8MF)G] = E[<p(F)H a (F, G)}, 

N 

Hi(F,G) = -Y,{G(DT F ,DF3) LHo>T) +T F (DG,DFJ) L2(0>T) -T F GLFi}, Mi G [l,iV], 

3=1 

H a (F,G) =if( aij ... )Q!m )(F,G t ) = H am (F,H^ air .. >am _^(F,G)). 

Also, for all q > 1, and all multi-index a, there exists (C,qo,qi,q2,ri,r2) only depending on 
(q,a) s.t. 

\\H a (F,G)\\ q < C||r F || go ||G|| gi>ri ||F|| g2 , r2 . (3.3) 

For the first part of the proposition we refer to Section V-9 of [21] , Concerning equation fj3.3[) . 
it can be directly derived from the Meyer inequalities on 1 1 Li* 1 || g and the explicit definition of 
H, see also Proposition 2.4 in Bally and Talay [3]. 

A crucial consequence of the integration by parts formula is the following representation 
for the density. 

Corollary 3.4 (Expression of the density and upper bound) Let F = (F%, ■ ■ ■ ,Fn) G 
(D°°) satisfy the nondegeneracy condition [ND]. The random vector F admits a density on 
M. N . Fix y G R . Introduce \/(u,v) G M 2 ,(/3q(u) = Iy> u , <4>\{v) = I v < u . For all multi-index 
j3 = (/3i, • • • , j3ff) G {0, 1}^ the density writes: 

N N 

p F (y) =nl\^(F)H a (F,l)](-lf\, a = (l,---,N), |0| := J>. (3.4) 

i=l j=i 

As a consequence of (|3.4p and (|3.3[) we get for all multi-index j3 G {0, 1} N : 

N 

3C>0, My)<Cn^(^r (i) ||H Q (^l)|| 2 , 7 « = 2- (m) . (3.5) 

i=l 

Proof Let B := U^Oi,^], Vi G [1, JV|, a t < h. Denote for all u G M, J («) := 
(-oo,u), Ii{u) : = [u,oo). Set finally, for all multi-index (3 G {0,1}^, Vy G R* *£(y) = 
A-rjv 7 / .•) Is(x)dx. Proposition 13.31 applied with a = (1, • • • ,N) and ^^ yields 

E[a a ^(F)] = E[^ B (F)H a (F, 1)]. (3.6) 
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Now, the r.h.s. of equation (J3T6J) writes 

N 
E[^ B (F)H a (F,l)} = E[ I B (y)dyH a (F,l)}= nl\l ytel0m) H a (F,l)}dy 

Jn&iifiiiFi) j b t=i 

N 

= n\{^%{Fr)H a (F,l)}dy. (3.7) 

Jb i=i 

The application of Fubini's theorem for the last but one equality is justified thanks to the 
integrability condition (|3.3p of Proposition 13.31 On the other hand, the l.h.s. in (|3.6p writes 

N 

E[d a * B (F)} = E[Hl FMaM (-lf>] = (-l)M / p F (y)dy. (3.8) 

Equation (J3.4J) is now a direct consequence of (|3,6p , ()3.7p , (|3.8p . Equation (|3.5p is then simply 
derived applying iteratively the Cauchy-Schwarz inequality. □ 

4 Malliavin Calculus to Derive Upper and Diagonal Bounds 
in our Examples 

4.1 Strategy and usual Brownian controls 

We here concentrate on the particular case of the process f)1.3|) (indeed the estimates con- 
cerning (|1.4p can be derived in a similar way). Since condition [H] is satisfied, assumption 
[ND] is fullfilled. It then follows from Theorem 2.3.2 in [32J that the process (X s ) s >q admits 
a smooth density p(t,x, .) at time t > 0. Our goal is to derive quantitative estimates on this 
density, emphasizing as well that we have different regimes in function of the starting/final 
points. 

To do that, we condition w.r.t. to the non-degenerate Brownian component for which we 
explicitly know the density. For all (£,x,£) G M + * x (M n+1 ) 2 we have: 

p(t,X,£) =p X l,n(t,Xl tn ,^ n )p X n + l(t,X n+ l,£ n+ l\X ' n = Xl, n ,X t ' n = £l >n ), 

( + c \ 1 ( I6-." ~ Xl >* 

P X l, n {t,X hn ,S,l,n) = JZ Tw5 eX P 



(27rt)»/ 2 * \ It 

We then focus on the conditional density which agrees with the one of a smooth functional, 
in the Malliavin sense, of the Brownian bridge. Precisely: 

p xn +l(t,X n+ i,£, n+1 \X ' =Xi tn ,X t ' = £i tn ) ;= p Yt {£ n+ \ - X n+ l), 



Yf.-- 



XlJ-^+Zl,nJ+Wy k 



du, (4.1) 



where (W u ' ) U £[o,t] i s f^e standard n-dimensional Brownian bridge on the interval [0, t\. The 
estimation of py t is the core of the probabilistic part of the current work. 
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We recall, see e.g. |34j . two ways to realize the standard n-dimensional Brownian bridge 
from a standard Brownian motion of M. n . Namely, if (Wt)t>o denotes a standard n-dimensional 
Brownian motion then 

(W u - jW t ) uem (1 = } (W^) uem , (4.2) 

Ut-u)f U ^A {l =\w»>% &m . (4.3) 



To recover the framework of Section I3.2| in order to deal with functionals of the Brownian 
increments, it is easier to consider the realization of the Brownian bridge given by (|4.3p . 



Remark 4.1 The process (W u ) ue \o jt ] := (Wt- U — Wt)ue[o,t] ^ a Brownian motion. More- 
over, the processes (W u — jWt) u e[o,t] an d ((* ~~ u ) Jo" iz~f)ue[o,t] are standard n-dimensional 
Brownian bridges on [0, t] , as well. 

For the sake of completeness, we recall some well known results concerning the Brownian 
motion and Brownian bridge. 

Proposition 4.1 Let q > 1, and (Wt)t>o be a standard n-dimensional Brownian motion. 
Then, there exists C := C(q,n) > s.t. for all t > 0, 

E[|WiH < C(q,n)t q/2 , E[ sup \W s \ q ] < C(q,n)t q/2 , 

se[o,t] 

E[ sup \W^\ q ] < C(q,n){t - r) q / 2 , < r < t. 
se[r,t] 

Moreover, there exists c := c{n) > 1, s.t. for all C > 0, and < r < t, 

/ I/-I2 

sup \W^\ > C] <2exp' h 



•e[T,t] V c(")(*-^) 

Proof. The first inequality is a simple consequence of the Brownian scaling. The second one 
can be derived from convexity inequalities and Levy's identity that we now recall (see e.g. 
Chapter 6 in |34j ) . Let (Bt)t>o be a standard scalar Brownian motion. Then: 

sup B u {1 = ] |£ s |,Vs>0. (4.4) 

ue[o,s] 

The third inequality follows from the first two and the representation (|4.2p . Eventually, the 
deviation estimates follow from (|4.4p as well. These deviations estimates can also be seen as 
special cases of Bernstein's inequality, see e.g. [34J p. 153. □ 

4.2 Some preliminary estimates on the Malliavin derivative and covariance 
matrix 

We now give the expressions of the Malliavin derivative and covariance matrix of the scalar 
random variable Yj defined in (|4.1|) and some associated controls. 
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Lemma 4.2 (Malliavin Derivative and some associated bounds) Let us set 

m(u,t,xi t n,£i >n ) := xiprj*- + £i,nf- Rewrite 



t ' *= L ,n t 



rO.tlk 



Y t = du\m(u,t,xi )n ,t,i yn ) + W® 

Jo 

du{\m(u,t,x ltn ,^ n )\ 2 + \W^\ 2 + 2(m(u,t,x lin ,Ci, n ), W^)} k ' 2 
J2 C l/2 I du\m(u, t, x hn ,^ n )\ k - 2i {\W^\ 2 + 2(m(u, t, si,„, &,„), W*>'>} 



i=o Jo 



(4.5) 



Considering the realization (|4.3p of the Brownian bridge, the Malliavin derivative of Yt (seen 
as a column vector) and the "covariance" matrix (that is in our case a scalar) write for all 
s£ [0,t]: 

fc/2 t 

D s Y t = Y. C k/2 ^|m( U ,t,x lin ,ei,n)| fc ~ 2i i{|VF u 'T + 2(m(n,t,x 1 , n ,ei,n),VK n ' t )r- 1 

4=1 ^ S 

fc/2 
t — U 



x 2-—— (W°'* + m(u, t, xi,„, &,„)) := J^ M *( S ' t, si,„, £i,„) 



t-s 

1=1 



ft 



'0 

Introduce now 



[ ds\D s Y t \ 2 . (4.6) 

Jo 



rt t — u 

M 1 (s,t,xi tn ,^ n ) := k du\m(u,t,xi, n ,£,i, n )\ k ~ 2 -—m(u,t,xi, n ,{, ltn ) 

J s IS 

+M/ ? (s,t,x 1>n ,a,„) := (Aff + M 1 fl )(a J t,zi,»,£l,n), (4.7) 

fc/2 

■R(s,Mi, n ,£i,n) := M ] R (g, t, xi, n , Cl,n) + y~]Mj(s, t, x 1>n , g 1|Tt ), 



i=2 



(4.8) 



Tit = f ds\(Mf + R)(s,t,x ltn ,^ n )\ 2 . 
Jo 

Set for allr € [0,£], 

Mr,* := f ds\M?(s,t,t ltn ,t hn )\ 2 , Mf.= M 0t t, 

K T>t := J ds\R(s,t,^ n ,Ci,n)\ 2 , K t :=K , t . (4.9) 

There exists C := C(k,n) > 1 s.i. /or a// r £ [0,i]: 
C-i(< - rj^lxi^l 2 ^" 1 ) + |6.«| 2(fc_1) ) < -Mr,* < C(t - r) 3 (\x 1>n \ 2 ( k -V + |e 1 , n | 2 ( fc - 1 )).(4.10) 
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Also, for all q>l, there exists C(k,n,q) s.t. 

n\Kr,t\ q ] 1/q < C(k,n,q)(t - rf(\x hn \V l^n]) 2 ^ 

x (t-r) f jt-r) 1 ' 2 \ 2(fc - 2) 



(4.11) 



Vk > 0, F[n T , t > KM T , t ] < c(n,k)exp (-K 2 ^ n \\j^^ ) , ( 1,12) 

for some constant c(n, k) > 1. 

Remark 4.3 From (|4.10p and (|4.1ip . it follows that 

2(fc-2) 



E 



M ./,<gw<:( 1+ i) ""'« 



w/ien |xi,n| V|£i )n | > Kt 1 ' 2 . For K := K(k,n,q) large enough, then the term Ait (correspond- 
ing to the Malliavin covariance matrix of a Gaussian contribution) dominates the remainder. 
This intuitively explains the Gaussian regime appearing in ii) of Theorem \2.1\ 



Proof. Assertion (|4,6p directly follows from the chain rule (see e.g. Proposition 1.2.3 in |32J) 
and the identity D^f = I s <ujz% , V(u,s) G [0,t] 2 deriving from (Ol) . 

Concerning (|4.10|) . we only prove the claim for r = for notational simplicity. Usual 
computations involving convexity inequalities yield that there exists C := C(k,n) > 1 s.t. 

Mt < Ct 3 (\x t , n \ 2 ( k -V + |ei,n| 2( ^ 1} ). (4.13) 

On the other hand to prove that a lower bound at the same ordre also holds for Mt one has 
to be a little more careful. 

W.l.o.g. we can assume that |£i,n| > |^i,n|- Indeed, because of the symmetry of the 
Brownian Bridge and its reversibility in time (see Remark 14. ip . if |£i n | < \ x i,n\ we can 
perform the computations w.r.t. to the Brownian bridge (W u ' ) u e[0,t] := 0^t-u)ue[0,t] using 
the sensitivity w.r.t. to the Brownian motion (W u ) u£ [ j] '■= (Wt-u — Wt)«e[o,tl- Note that 
|6,n| > kl,n| => |6,n|oo > ^ |xi,n|oo- Let i Q G [l,n] be the index s.t. |£i,n|oo : = |& |> tlien 
|Ciol — ^72l x iol- Let us now write 

,2 /"* , / /"* , i i » Mfc-2 ft - u u.\t — u\ 

M t >k J as I du\m{u,t,xi in ,^i tn )\ I — — x io + -& I t^— I • 

Observe now that for s > ^" 2 t we have that \/u G [s,t], ^j^£j + f £*o ^ as the s iS n °f Cio- 



n 

Hence, 



M t >k 2 I , ds[ I du 



t — U u 

7 x io ~r Tsio 
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fe-i , N 

t — u 



(4.14) 



,1/2 



Now, for s > t ( i/" , 2 _i ) , we have for all u S [s, t]: 



t — u u 

7 x io i Ts«o 



k-l 



> ( - 



> \ii 



k - l \ii \ k ' x ft 



2 k-2 



U 



|fc-l 



n 



1/2 



t 

k-l 



k-l 



ik-1 



J »ol 



2nV2 + i 



(4.15) 



Equation (|4TTTJj) thus follows from (|4TT4"j) . (|4TT5l) and ([4TT3J) . 

Concerning the remainders we get that there exists C3 := 63(77, k),C4 := C^{n,k) s.t.: 

|Mf( S ,t,X ljn ,6,n)| 2 < ^(t-s) 2 !^!^" 13 SUp K'f I6.nl" 2 , 

u<=[s,t] 

sup |<>f (2i " 1) |6,n|- 2(2i - 1) 

u£[s,i] 



+ SUp I^T^I^nl- 2 ^)- 
ue\s,i\ J 



(4.16) 



From (|4.8|) and a convexity inequality, we derive \R(s,t, xi :n , S,\ :n )\ 2 < 2{\Mf-(s,t,xi !n ,S,i :n )\ 2 + 



Yli=2 \^i{ s ^,xi, n ,^i jn )\ 2 ). Thus, from (|4.16p . (|4.9|) . we obtain that there exists C(k, 
s.t. for all r € [0,t\: 



n. 



. fc/2 



E 



[\n T , t \^ < c(k,n, q )(t-Tf\^ n \ 2 ^\^n\ sup K-f^- 1 )!^.^- 2 ^- 1 ^] 179 

l i=1 ue[r,t} 



fe/2 



+ £l[| SUp IVeT^kl.r 



-2(i-l)|gil/g 



_2 u£[T,t] 

which, thanks to Proposition 14.11 gives (|4.11|) . 
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On the other hand, from (|4.10p and the previous convexity inequality for R we get: 

n^M T ,t < Kr,t\ < P C-\t-Tf\^ n \ 2{k - l) K<^ f ds\Mf(s,t,x hn ,^ n )\ 2 

fe/2 ,t -1-1 

+Y^ / ds\Mi(s,t,x ltn ,Ci,n)\ 2 

A — O J T -I - 

\ C -l (f _ r) 3| ein |2( fe -l) K <y d8 | A f^( a ,*,Xi,n,6,»)l 



fc/2 

+£ 

i=2 



) C-\t-Tf\^ n \^ k - l ^K< f ds\Mi(s,t,x lin ,^ 



HA6t 

< 



fc/2 

-E 

i=2 



l) C-Ht-rf^f^KK^t-rf^nf^ sup K'f|6,n|- 2 
K J 6 ue[r,t] 

o\ 2 ^-y ( / C11T > IW°'*I \ 2(2i-l) / c ™ IW '*! \ 2(i-l) 

2 \ x C 4 J /sup ue [ r4 ] \Wu \\ v v / sup„ s [ T)t ] I kv u |\ v > 



Equation (|4.12|) then follows from Proposition 14.11 □ 

4.3 Control of the weights 

Now to exploit Corollary [33] to give estimates on py t we need to have bounds on the Malliavin 
weights. Formula (|3.4p involves two kinds of terms: the inverse of the Malliavin Matrix and 
the Ornstein-Uhlenbeck operator. Lemma [4. 21 provides tools to analyze the Malliavin matrix. 
Concerning the Ornstein-Uhlenbeck operator we will rely on the chaos expansion techniques 
introduced in Section [3JJJ 



4.3.1 "Gaussian" regime 

In this section we assume that |xi n | V |£i,n| > Kt 1 ' 2 , for K := K(n,d) sufficiently large. 
That is we suppose that the starting or the final point of the non-degenerate component has 
greater norm than the characteristic time-scale t 1 ' 2 . In this case, we show below that the 
dominating term in the Malliavin derivative is the one associated to the non-random term 
Mj° in (|4.7p . This term corresponds to the Malliavin derivative of a Gaussian process. This 
justifies the terminology "Gaussian" regime. 

In order to give precise asymptotics on the density of Y t , the crucial step consists in 
controlling the norm of Ty t '■= 7y in L q {Q), q £ [1, +oo) spaces. 

Lemma 4.4 (Estimates on the Malliavin covariance) Assume that |xi >n |V|£i )n | > Kt 1 / 2 . 
Then, for all q G [1, +oo) there exists C„0~J]:= C ' w^(n,k,K) > 1 s.t. 



C g ™ 



(| s |2(fc-l) + |£ |2(fc-l) )t 3 



*<PYX< 



c. 



( | Xln |2(fc-l) + | ei>n |2(fc-l) )t 3- 
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Proof. As in Lemma 14.21 we assume, without loss of generality, that |£i«| > |#i,n|- To give 
the L q estimates of the Malliavin derivative we recall the definition of Ait given in (|4.9p . and 
we use the following partition: 



n\ry t \ q ] = E E [l r ^l ?I r y6 [^ ffetili 



< 



M, 



m>l 



M 4 



4m J 



(4.17) 



Equation ()4.10p in Lemma 14.21 provides us with an useful bound for Ait- We next give 
estimates of P \yy t < ^1 , m > 1 in the spirit of Bally pQ. 

Introduce t m := inf{t> E [0,t] : A^«,i < Ait/m}. We first show that there exists mo 6 N 
and C := C(n, A;) such that t m > i(l — Cm~ 1 ' 3 ) for all m > mo- From (j4.10p we obtain 

t m > mi{v€[0,t]:M v ,t<Ct 3 (\x 1 , n \^ k -^ + \^ n \ 2(k -^)/m} 
> M{v G [0,t] : Mv,t < 2Ct 3 |6,n| 2(fc " 1) M} =■ ?m, 

recalling we have assumed |£i, n | > |^i,n| f° r the l as t inequality. Equations (|4.14p and (|4.15|) 

1/2 

also yield that there exists C2 '■= C^in, k) s.t. for all v > ^ , t, 

^.t^Ca^-t;) 3 ^! 2 ^- 1 ). 
► £ as m — )• +00, then there exists m such that t m > 



Note that L 



71 



1/2 



w i/2 +2 -i * for eVel T 

m > m, and the above inequality holds for every v G [t m ,t]. Set C := (2C/C2) 1 ' 3 , and 



1-^3, 



^o = \C J V m. For every m > mo we have that: 

rt ka 

< P / ds\(MP + R)(s,t,x 1<n ,Z lin )\ 2 <— - 
Ut m 4m. 



~ 4m 



< 



< 



ds\Mi '(s>t, x 1>n , Ci, n )| 2 - / c?s|i?(s,t,j;i iri ,^i )n )| 2 



< 



X 



1 771 it 



< Tit 



< c(n, &) exp 



21/3 



|6,n| m 
16Cc(n, fe)t i ' 



4m 

(4.18) 



using (|4.12p for the last inequality. Plugging this control into (|4.17j) , using once again (|4,10p 
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we derive that there exists C3 := 63(72, /c), (64,65) := (64, C§){n, k, q) s.t.: 



|gi,n| r n 
t 



2™ 1/3 



, . ^ „y 4 



t 3 l&,»P- 1 V Vlft,~l (2( *- 1)+6 » 



C S (8C)« 



E 

m>mo 



V3IA _|2 



m V<»|£ 



3g 



l.n 



exp -C 3 x 



|Ci,n| m 

t 



2 1/3 



I „."*j„Y + 



C, 



t 3 !^! 2 ^- 1 ); i£i,„p +4 >* 



<a 



+ 



^|£ |2,(fc-l) |^ lin |(2fc+4), | 6>ri | 



< 



£ exp -C, 

m>mo \ 

c 5 



-iiei,n|W/3^ 



^ |2,(fc-l) 



1 + 



(3(9+1) 

|£l,n| 6 <* +1 > 



which for |£i n | > Kt 1 ' 2 gives the upper bound of the lemma. 
Let us now turn to the lower bound for ||r>^ ||z,p(]p) - Write: 
E\T q Yt ] > E[T q Y I^ Yt < 3Mt ] > 77r L x -^[ lYt < 3M t ] > j^j^i I F'h, - :]>f/| ). 



(3M 



(3M 



From equations gSJ-gS]) one has Ppy^ > 3M] < P[2_M t + 2K t > SM t ] = F[Ht > |M]- 
Now, from Lemma IO equation (|4.12|> . one gets P[7y t > 3M t ] < c(n, k) exp (— jJ^L ) . 

Therefore, for |£i, n | > Kt 1 ' 2 and if large enough, we get E[Ty ] > w^n?) which thanks to 
(j4.10p completes the proof. □ 



Controls of the weight for the integration by parts. 

From Proposition 13.31 and Corollary 13.41 we derive 

PY t (tn+l ~ Xn+l) = E[H t I Yt >t n+1 -x n+1 ], 

H t = -(DT Yt ,DY t ) L2m + T Yt LY t = 7y 2 (D lY , DY t ) L 2 m + T Yt LY t 

■= Hl + H 2 t , (4.19) 

using the chain rule for the last but one identity. 

We have the following L 9 (P), q > 1, bounds for the random variable Hf. 

Proposition 4.2 (Estimates for the Malliavin weight) Assume that |xi )n |V|£i )Tl | > Kt 1 ' 2 
for K large enough. Then, for all q G [l,+oo) there exists C „ra~2b= ^ w^n,k,K) > 1 s.t. 



Ht\\ q < 



c aWm 



(Knl^-D + I^nl^- 1 ))* 3 / 2 ' 
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Proof: Control of H\ . From (|4.19p we get for all given q > 1, 

\\Ht\U ■■= E[7y, 2 *l(O7n.O>'.>L«(0,,)n 1/ " < Eb^FMl (°TC , W.)x»(o,.)| 2 «] 1/2 ' 

- -(fc.Ai^ ' i " > » i,jy '' gMn "' 

- ^(l6.l^'° 1 ,„l^--»)^ l|DTCl ''^ 1 ' /4 ° E[|OY ' l ''^ 1 ' /4g - (4 ' 20) 

using Lemma 14.41 for the last but one inequality. Now, from equations (|4.6p . (|4.8j) . using the 
notations of Lemma 14.41 

< (2 4 *- 1 [M 2 t q + E[TZ 2q ]}) 1/4Q < 2 1 - 1 ^ [m] 12 + EiK?] 1 ^} . 

On the one hand equation (|4.10p in Lemma 14.21 readily gives M t < Ct 3 ' 2 (\xi >n \ + 
|^i,n| fc_1 )- On the other hand, equation (|4.12|) of the same Lemma yields 



Hence, there exists C\ := C\(n,k,q, K) s.t. 

n\ DY t\% m ] 1/4q = E h? t ] 1/iq < C 1 t i l 2 {\^ n \ k - l + \x^ n \ k - 1 ). (4.21) 

In order to get a bound for ||flj || g , it remains to control E[|Z>yyJ L 'L 0f s] 1 / 4 ' ? . Equation (|4.6p 
and the chain rule yield that for all ui G [0, t], D U2 jY t = 2 J du\D U2 D Ul Y t x D Ul Y t . We get 



t ft 

2\4<jnl/8g 



E [l^y t lS (0jt) ] 1/49 < mi^mj d Ul J du 2 \D U2 , Ul Y t 

< Cs^da.nl^ 1 + Ixi^^EOZ) 2 ^^^] 1 /^ (4.22) 

C2 := C2(n,k,q,K) using (|4.21|) for the last inequality. 

With the notations of equations (|4.6p . (|4.8p we set for all u\ £ [0,t], 

fe/2 fc/2 

i=l i=l 

for simplicity. 

Observe now that for all i 6 [2, fe/2], «2 £ [0,£], 



D U2 M i {u 1 ,t) = C l k/2 / ^|m(M,x lin ,£i,„)| fc - 2i {|W^<f + 2(m(M,xi,n,6,n), W^ 4 )}* 

x2z ^"^ {2(i - 1){W^ + m(v, t, x 1>n , &,„)) ® «»' + ro(v, t, si,„, &,„)) 

(t-ui)(t-u 2 ) 

+ { I W°'f + 2(m(«, t, si, nj 6,n), <'*)} /«} , 

fc-2 (* " V) 2 



D U2 Mi(ui,t) = k / df|m(u,i,xi.„,^i.„)| — — -I n . 
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Prom the above equations, assuming once again w.l.o.g. |£i,n| > |^i,n|i the arguments used 
in Lemma 14.21 yield: 

nD 2 Y t \% m)2) } 1 ^ < CE[\ J 2 d Ul du 2 (t - u x V u 2 ) 2 |£i,n| 2(fc - 2) 

fe/2 

X (l + V|6, n | 4(1 "' ) sup i^.fCi-Dji^i/s, 

fe/2 

< Ct 2 |6,n| fc - 2 (1 + J] I6,n| 2(1 - J) 1E[ sup |w£.*|ie«P~i)]i 

t=2 ue [ '*] 



m ED 7 / /1/2 \ 



), 



where C := C(n,k,q) may change from line to line. Recalling that |£i n | V |xi,n| > Kt ' 2 we 
obtain 

E[|£> 2 r t |g ((0jt)a) ] 1 / 8 « < Ct^i,*!*- 2 , C == C(n,fe )9) if). 

Plugging the above equation into (|4.'22p we derive that 

n\Djp t \ L , m ] 1/4q < ^ 7/2 i6,ni 2fc - 3 , 

which together with (|4.2ip and (|4.20p . eventually yields 

Control of if 2 . From (|4.19p and Lemma I4.4I for all g > 1, we get 

KII, < E[|r y j 2 Y /29 E[|Ly 4 | 29 ] 1/29 < ^77^T^^^ E [l Ly 'l 29 ] 1/29 - ( 424 ) 

Now, since LYj = <5(-DYt), the idea is to provide a chaotic representation of DYf. To do 
that, we use Proposition 13.21 (Stroock's formula see [38]). For a given u\ G [0, £], recalling 
fc/2 

D Ul Yt := > Mi{u\,t) where Mj(ui, t) E R n is a random contribution involving Wiener chaos 

i=l 
up to order % — 1, one has: 

2t-l 
M 4 (m,i) = E^Cui^JlH- J^ J,(^j(., ui,t)), 

z=i 

Ii(gl(.,ut,t)) ■= ••• / 9i(vi,--- ,vi,u 1 ,t)®dW Vl ®---®dW vl , 

Jo Jo Jo 

gj( Vl ,--- ,v l: u 1: t) := E[D Vlt ... iVl M i (u 1 ,t)]G(R r T (l+1 \ (dW Vt ® ■ ■ ■ ® dW Vl ) € ((R n )® 1 )* . 
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fe-1 fe/2 

Hence, An^t := go(ui,t) + > /;(g;(.,ni,t)), where 5o(«l,*) := sMMj(ui,t)] and for all 

1=1 i=i 

fe/2 

Z e [l,fc-l], 5j(ui,--- ,uj,wi,t) := ^ ^K^i, • • • ,t>j,ui,i), so that 

i=[l/2]+l 

LY t = j g (ui,t) ® d^ Ul +X)/j(fli-i (.,<)) :=X)/,(gi-i(.,t)). (4.25) 

•^° «=2 J=l 

Similarly to the proof performed to control IE[|"7^t | 2<3r ] 1 / 4<? of (|4.12|) in Lemma 14.21 we obtain 
that there exists C := C(n, k) s.t. for all I £ [0, k — 1] and for all (vi, ■ ■ ■ ,vi,u\) S [0, t] . 

\gi( Vl ,--- , Vl , Ul ,t)\ < Ct(|^,n| Mm) + |xi, n | Mm) )- (4.26) 

Therefore, 

E [| L y t |2«]l/2, < c jy+J/2 ( |£ i>n |fc-J + | Xlin |*-I) 

Z=l 

where C := C(n,k,q) may change from line to line. Recalling that |£i n \ V \xi n \ > Kt 1 ' 2 , we 
derive from (|4.24p that there exists C2 := C2(n, fc, g, X) s.t. 

which together with (|4.23p and (|4.19|) completes the proof. □ 

4.3.2 Non Gaussian regime 

We now consider the case |xi n | V |£i,n| < Kt 1 ' 2 , which corresponds to a diagonal regime 
of the non-degenerate component w.r.t. the characteristic time scale. It turns out that 
the characteristic time-scale of the density PY t {in+\ — %n+l) is t 1+k ' 2 . Indeed, we have the 
following result. 

Proposition 4.3 (Estimates for the Malliavin weight in Non Gaussian regime ) Let 

K > be given and assume that \x\ n \ V |£i,n| < Kt 1 ' 2 . For every q > 1 there exists 
C aW3 ■= Qofa fc > K ) s± 

ilwil < Cq 



t l+k/2- 
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Proof. For t > write: 



Y t 



Thus: 



o 



t — u . u rT ..„ 
asi,«— — +£i,n~ + W£ 



o.t 



du 



a+fc/2 






iV2 



*V2 *l/2 



du 



-: t l+k ' 2 Y\ 



PYt (in + i-x n+1 ) := -d in+ J[Y t > Z n+ i-x n+1 ] = -d &i+1 P[Y{ > e " ' " ■''"-'- 



t l+k/2 



1 /£n+l — X n+ i 

>Pvt( — JUTT? — )■ 



From Corollary 13.41 (Malliavin representation of the densities), we obtain: 



PY t {t,n+l ~ Xn+l) = E[H(Y t ,l)I Yt> ^ n+1 - 



Zn + lJ 



E[H(Y\,l)I^^ n+1 - Xn+1 \ 



t l+k/2 



~Yl>- 



jl+fc/2 



' E[H(Y\,l)I Yt>in+1 _ Xn+1 ] 



t l+k/2 



so that H t := H(Y t , 1) = t~^ +k ^H{Y\, 1) := H 1+fe / 2 ).Hf a . Hence, for all ? > 1, 



U s' TJ 1 

1^*119 - ^TTfc72" l|jH i H«- 



(4.27) 



Now, as a consequence of the Brownian scaling we get (-t$-)ug{o,i] = (W«' ) u e[o,i] so 

I x l," I w I 51,! 



du. Recalling that |^| V |S^| < K we 



that P* (1 = } t 1+fc / 2 Jo 1 |ff (1 - «) + §£u + W. ' 1 

derive that the usual techniques used to prove the non degeneracy of the Malliavin covariance 

matrix under Hormander's condition (see e.g. Norris [31 J or Nualart [32 1) yield that there 

exists C q := C q (n, k, K) G R+* s.t. \\H{ 1 \\ q < C q which from (J4T2?)) concludes the proof. The 
crucial tool here is the global scaling. □ 



4.4 Deviation estimates 
4.4.1 Off-diagonal bounds 



From the Malliavin representation of the density given by (|4.19j) . to derive off-diagonal bounds 
on the density, it remains to give estimates on P[Yj > £ n +i — a?n+l]- 

Lemma 4.5 (Off-diagonal bounds) LetU k (x,^) : = £ n+ i - x n+ i - 2 -+ r (|xi, n | fc + |£i,n| fe )*, 
and assume that Ut(x,£) > 0. Then, there exists Cfe~%\ '-= C j^5] (n, k) s.t. 

(i) If \x\ >n \ V |£i jn | > Kt 1 / 2 for a given K > 0, 

m > e, +1 - x n+l] < to {ex P (-c ^ (| |fc 5|7^ B 2 | fc -i) 



) 2 £ 3 



+ exp 



1 |a>i,n| 2 + |f 



l,n 



fc/2 



E ex p Hfck 



u t k (x,0 lfi 



{\xi, n \ k ~ 2i + l^nl*- 2 *} 1 /** 1 ^/* 



24 



(ii) If |xi )Tl | V |^i,n| < Kt 1 ' 2 for the same previous K , 

m > * n+1 - x n+ i] < ^{exp (-^%|2) +i>p (-q yK' ) }• 

Proof. We only prove point (%), the second point can be derived in a similar way. According 
with (|4.5|) . we first decompose Y t as 

5t= / |m(u,t,xi, n ,Ci,n)| fc ^ + M t fc (a;i )n ,^i in ) + J R t fc (xi in ,6,n) 

JO 



where 



Mt(xi >n ,^i,n) ■= kj \m(u,t,xi tn ,Ci,n)\ k 2 {m(u, t, xi >n , Ci,n), Wu'^du, 
RHxi,n,Ci,n) := | £ |m(u,t,si >nj $i, n )|*- 2 |wS'*| 2 d« 

fe/2 t 

+E C fc/2 / |m(n,t,x 1 ,„,6,n)| fc - 24 (2(m( U ,t,x 1 , n ,6,n),^ i ) + W| 2 )*du, 

i=2 ^ 

then we have 

F[^t > £n+l " X n+l ] = F[M t fc (x ljn ,ei,n) + i? t fe (x ljn ,ei,„) > 
£n+i -a^n+i -J |m(u,t,xi jn ,^i in )| fc (iu]. 

Note that all the terms in R^ (xi >ri , £i,n) have characteristic time scales that are in small time 
negligible with respect to the one of the Gaussian contribution M t fc (xi jn ,^i jn ). Moreover 

M*(x!,m£l,n) < 2 fc - 2 (|x 1 , n | fc " 1 + l^l.nl^ 1 )* SU P \W^\ =: M t fc (x ljn , 6,n). 

us[0,t] 

Since by assumption ?7 t (x, £) < S, n +i — x n+ \ — L \m(u, t, xi,n>£l,n)| du, one gets: 

W[Y t > £ n+l - x n+1 ] < F[(M t k + i2*)(xi, n ,6,n) > ttffoO] < 

P[2M t fc (x 1 , n ,6,n)>C/ t fc (x,0] 
+ P[(M 4 fc + J R t fc )(x 1 , n ,6 : n) > ^,£)] 1/2 

xP[i? t fe (x ljn ,ei,„) > M t fc (xi, n ,ei,n)] 1/2 - (4.28) 

Standard computations, similar to the ones performed to prove the deviation estimate in 
(|4.12p in Lemma 1431 give that there exist C\ := C\(k), C 2 := C 2 (n, k) > 1 s.t. 

F[RUxi, n ,Ci,n) > M t k (x 1>n ,£ 1>n )} < (k - 1)P[ sup |W°-'| > Ci{|xi,„| + |£l,„|}] 

ue[o,t] 

|2 , It |2" 



< C 2 exp I -C 2 1 

(kLnl^ + l^nl*- 1 )* 3 



P[M t fc (x 1 ,„,ei,n)>^(x,0/2] < Csexpf-^ 1 - -^l|^— — ) . (4.29) 
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On the other hand, we have: 



[(Mf + fc.U^^^)] 



Now, 



< P[M 4 fc (x!, n ,6,n) > lu t k (x,0] +F[^ fc (x lin ,6,n) > lu t k (x,$}. (4.30) 



P[fifK^l,n)>^f(^)] 
o(fc-3)V0k i 

<n 9(h n {|si,n| fc - 2 + lft.nl*-- 2 }* SUp \W^\ 2 >——-U t k (x,0} 

tk-i-2 



ofc+i-2 1 ^i 



ue[o,t] 

fe/2 

:=P 1 + ]T(Pj+Pj). 



i=2 



(4.31) 



From Proposition 14.11 one gets that there exists C3 := Cs{k,n) > 1 s.t. 



{|zi,„l*- 2 ' + I6,,.l t - 2 '} 1/ '* 1+1/ V ' 

H £ ^^(-^' {K^rlfe^ )- (4 ' 32) 

Hence, plugging (jQ2j) in (OB we derive the claim from (JOB . QUI) . (109]) and (jOHD-D 

4.4.2 Auxiliary deviation estimates 

Still from the Malliavin representation of the density given by (|4.19p . when £ n +i — x n+ \ is 
small, that is when for the degenerate component the starting and final points are close, we 
have to give estimates on F[Y t < £ n +i — x n+ i] (small and moderate deviations). 

Proposition 4.4 There exist constants (01,02) := (ci,C2)(n,k) s.t. for all (xi n ,£i n ) £ 
(R"\{0}) 2 , en+l > x n+x and t > 2 fc + 3 ,/jl^V ; 

/ I |2+fe , It |2+fc\ 

m < £n+i - x n+1 ] < Cl exp -c 2 1 -^ -tlSLn! . (4.33) 

V 4n+l - X n+ i J 
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For a given K > 0, if t > [{^ n+1 - x n+1 )\{UK) k ] 2/{k+2) , and \xi, n \ V |£i, n | < Kt 1 / 2 , then 
there exist (ci,c 2 ) := (c"i,c 2 )(n, k,K): 



W[Y t < £n+i - x n+ i] < c\ exp -c 2 



^1+2/fc 



(Cn+l - £n+l) 2/fe / ' 



(4.34) 



Proof. We first begin with the proof of (|4.33[> . As in the previous sections, we can assume 



w.l.o.g. that |xi, n | > |£i,n|-For s G [0, £], we define X s := xi jU ^- + £i,n§ + WV (where 
(W s ' ) se [o,t] is a standard n-dimensional Brownian Bridge on [0,£]), so that Y t = f \X s \ k ds. 
Let us also set ri^ i/ 2 := inf{s > : |X S | < |xi in |/2}. Consider now the event A := 
i T \x ln \/2 < 2 k \x~ X \k +1 } an d denote by .A c its complementary. Observe that F[f Q \X s \ k ds < 

e„+i - x n+1 ,A c ] = P[/ ^ (%^) ds < Jj |X s | fc ds < £ n+1 - x n+1 ,A c ] = 0. Thus, 
t < U+i - Xn+i] = V[Y t < U+i ~ x n+1 ,A] < F[A]. Now 
F[A] < P[ inf |X S | < |a?i,„|/2] 

t- s 



se[o,2 k in + 1 x " +1 i 



< P[ inf 



2-1, n , ~r Sl,n , 



+ inf (- W|) < \xt, n \/2] 

se[o,2 k in + 1 x ™ +1 ] 



< P[|xi, n |/2+ inf (- -){|xi, n | + |ei,„|} 

sgp^ ^l 1 Xn k +1 ] l 



sup iw;- 4 ! < o] 



»e[o,2* g "+ 1 "" +1 ] 
ki, ra r 



< P[|x lin |(l/2 



2 k+1 ^ n+1 -x n+1 ) 



\Xl n » 



< 



se[o,2 



sup |w fl 0i '| 

fc Sn + l-^Ti+l i 
Ui ! fe J 



< P[|xiJ/4< sup 



IW. 



O.ti 



s6[0,2 ft 



Jn + l~ x n+l l 



recalling Ixi n | > |£i J and t > 2 fc+3 ^"| 1 x ,1 +1 for the last two inequalities. From Proposition 
14.11 we obtain: 



[Y t < £„+i - x n+1 ] < F[A] < ci exp -c 2 



Fi, 



|2+Ai 



'Cn+1 — %n+l, 
which from the assumption \xx n \ > |£i,n| gives (|4.33|) up to a modification of c 2 . 

Let us now turn to (|4.34j) . Introduce i^g(i) := j IW ik<fl(£ -^ \ds for a parameter 
/3 > to be fixed later on. Define the set Ap := {I pit) > t/A}. Observe that 

ft 



[ \X s \ k ds <£ n +i-x n+1 ,A 



[/ ^.I^K^-^ol^-l*^^ / \X s \ k ds<U + i-x n+1 ,A c ^ 



< F[P{Zn+i ~ x„+i)3t/4 < / |X s | fc (is < e„+i - x n+1 ,A C s\. 

Jo 
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Choosing (3 = J^ we get from the above inequality P[J" |X s | fc ds < £, n +i — x n+ i,A^] = 0. 
Hence, 



ft 
< 



/ \X s \ k ds<ti n+1 -x n+1 ]=F[[ \X s \ k ds<^ n+1 -x n+1 ,A ± ]<F[A ± ] 
Jo Jo 3t 3t 

n [ U , k< 3j, n+1 - Xn+l) ds > t/4] < P[ / 1^ s(in+1 - Xn+l) ds > t/4] 

JO * s| - 4t Jo ' sl - 3i 

^+B°A<c(^t,k) dS > t/4] ' C ( X '^*' fc ) := ^ 4^ J ' 

< P[^ Vi^+^f +i#*|<c(*.CA*) <to > ^ 8 1 + P t/ /2 W+6f +**|<c(^A*) <to > */ 8 l 
:=Pi+P 2 , (4.35) 

where (P s ' ) s e[ot] stands for a one-dimensional Brownian bridge on [0,i\. Observing that 
(B s ' ) := (B t L s ) S £[o,t] is also a Brownian bridge, we get that 

rf/2 

= p[y o ^i| Xlf+ft *- + B0.*i< C (^, fe) ^>V8]- 

Since we assumed |xi| V |£i| < Kt 1 ' 2 , \x±\ and |£i| have at most the same magnitude so that 
Pi and P2 can be handled exactly in the same way. Let us deal with P\. The occupation 
time formula for semimartingales (see Chapter 6 in |34j ) yields 

t/2 r-c(x,i,t,k) 



where LL stands for the local time at level z and time t/2 of the process (xi^p + £i| + 
B s ' )s6[ot]- From the definition of Pi in (J4.35P : 

Pi < P[ sup L z t/2 x2c(x,^t,k)>l] 

ze[-c(x,t,t,k),c(x,£,t,k)] ° 

= P[ sup Ii /2 > * - ], (4.36) 

2fc ^ ^72 • ^172 1 

where L 1 / 2 stands for the local time at level z and time 1/2 for the scalar process 

(y \ ( Xl ft \ . & 1 B °A (law) f Xl /1 \ 1 & 1 rO,A 

(*«)«e[o,l] : = Ui7^( 1 -«) + ^«+^72 = (^( 1 - u ) + ^ n + jB « ) 

V / u6[0,l] / US [0,1] 

The last equality in (|4.36p is a consequence of the scaling properties of the local time. From 
Tanaka's formula for semimartingales L^m = |-^i/2 ~ A~ l-^o ~ A ~ fo sgn(X s — z)dX s . 
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R ' n 1 

Denoting with a slight abuse of notation (-fffi) U £[o t i] = {B u ' ) ug [o,i], we have the following 
differential dynamics for X u : 



dX u = -j- — du + dB u ' = -^ du + dB u , 



tV2 



l-U 



where (-Bn)ne[o,il i s a standard scalar Brownian motion. 

Therefore, from equation (|4.36|) and the usual differential dynamics for the Brownian 
bridge: 



P < Tur l^l ~ x ll i ir0,1| 



+ 



ze 



sup 

c(x,£,t,k) c(x,£,£,fc) -i 
?72~ ' t l/2 l 



1 'sgniXs-zX-^-^ds + dB*' 1 : 



t 1 / 2 



> 



t 



1/2 



8c(x,£, t, k)' 



< 



■16 -xi 



1/2 + l-°l/2l + 



i 1 / 



sup 

g,t,fc) I 

tV2 > t l/2 



sgn(X s -z)di3 a | > 



t l/2 



1/2 |r0,1| rl/2 

LfcO i kj 1,1 1/ I J UWli l J Jk B A^ * U/JJ g I ■- . , 

i! 1-s . C (x,g,t,fc) c( a ,e,t,fc) 1 Jo 16c(x, £, t, k) ' 



<P[2ET + 3 sup \B^' l \ + sup 

*6[- 



1/2 



sgn(X s -z)dB s | > 



t l/2 



se[o,i/2] - r r gfeljM) gfelAM i -/o 16c(x, £,*>&)" 



1/2 ' t l/2 



Now from the definition of c(x,£,t,k) in (|Q5]I . for t > [(£ n+ i - x x+ i)l(64K) k ] A + ' one 



fc1 2/(fc+2) 



has 



t l/2 



16c(x,g,t,fc) 



2K > nn fJ^,.s . Thus 



32c(x,£,t ] fc) ' 



t l/2 



+1 



sup 

|.t,fc) < 

*l/2 ' tV 2 



Pi<P[3 sup \B^\> ^ f 

»6[0,i/2] 64e{x,£,t,k) 

1/2 _ f l/2 

sgn(X s - z)dB s \ > 



. c r c(x,g,t,;.: .,,;,.; ./(i 64c(x,6<,fc) 



Setting for all t G [0,1/2], M t := J * sgn(X s - z)dB s , M t := B {M)t = B t (i.e. B is the 
Dambis-Dubbins-Schwarz Brownian motion associated to M). Hence, from Proposition 14.11 
we derive the announced bound for P\. Since P2 can be handled in a similar way, the claim 
then follows from equation (|4.35p . 

4.5 Final derivation of the upper-bounds in the various regimes 

In this section we put together our previous estimates in order to derive the upper bounds of 
Theorem 12. II in the various regimes. 



29 



4.5.1 Derivation of the Gaussian upper bounds 

In this paragraph we assume |xi,n| V |£i, n | > Kt 1 ' 2 for K large enough. We also suppose 
l€ " + ^7 2 !g ( j' g + 1 1 "jf- ( i+jgf)i-J" |fc)l < C where c := c(k) = 2+ ^ and C is fixed. From CorollaryE! 
(representation of the density) , Proposition 14.21 (controls of the weight in the integration by 
part) and Lemma 14.51 (deviation bounds), we have that there exists C := C(n,k,K,C) > 1, 
s.t. setting U^{x,^) := S, n +i — %n+i — ^rr(|^i,n| fe + |£i,n| fc )£ as in Lemma [431 one has: 

U ^ S 06XP l 2t _ (|xi,n| fc - 1 + |6,n| fc - 1 ) 2 tV ,, „-v 

**>*'*> ~ W 2+ 3/ 2( | Xl ^|,-l + | 6 ^| fc -l) • ^ 37 ) 

Remark 4.6 The above result means that the Gaussian regime holds if the final point £i >n 
of the degenerate component has the same order as the "mean" transport term mj(x,£) := 
x n+i + x+rd 3 ' 1 .™! + l£i,n| )t (moderate deviations). A similar lower bound holds true, see 
LemmaV 



4.5.2 Derivation of the heavy-tailed upper bounds 

We here assume ^ n+ Xi2n +1 ~^ihf + }k\^ > C where c := c(k) = 2 + ^ and C is as in 
t d/ ''(|a;i,nr -1 +|ci,nr _1 ) ~~ fe+1 

the previous paragraph. 

If \xi }n \ V |Ci,n| < Kt 1 ' 2 (K being as in the previous paragraph), then Corollary 13.41 
Proposition 14.31 and Lemma 14.51 yield that there exists C := C(n,k) > 1 s.t. 

P(t,x,0< t(n+fc)/2+1 exp I C tl+2/k I. (4.38) 

On the other hand if \xi >n \ V |£i, n | > Kt 1 ' 2 , then Corollary 13, 4( Proposition 14.21 and Lemma 
14.51 yield that there exists C := C(n, k) > 1 s.t. 

,, ^ , C / lfr,n-Si, n | 2 X-l UtHx,0 2/k \ 

W> X >U - t n/2+3/2^ Xln \k-l + \£ hn \k-l) 6XP ^ 2* tl+!»A J 

- K k-l t (n+k)/2+l ^ I 2t t X + 2 / k ) ' 

Hence, up to a modification of C, the control given by (|4.38p holds for all off-diagonal cases. 

4.5.3 Moderate deviations of the degenerate component 

In this paragraph we suppose < £ n +i — x n+i < Kt 1+ ' 2 , for K sufficiently small. This 
means that the deviation of the degenerate component is small w.r.t. its characteristic time 
scale. From Corollary 13.41 Propositions 14.21 and 14.31 and Proposition 14.41 we derive similarly 
to the previous paragraph that there exists C := C(n,k,K) s.t. 

r PYH ( Ifr."-*!."! 2 _ r-l f kl.n| 2 + fc + l€l.n| 2 + fc , f! + 2 / fc 1\ 

*..,« < H " 21 jj^r «— )a " )J .(4.3 9) 
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4.6 Gaussian lower bound on the compact sets of the metric 

We conclude this section with a proof of a lower bound for the density on the compact sets 
of the metric associated to the Gaussian regime in Theorem 12.11 A similar feature already 
appears in the appendix of [TU] . 

Lemma 4.7 Assume that \xi n \ V |£in| > Kt 1 ' 2 , K > Kq := Ko(n,k) and that for a given 

C > we have ^"V^* 1 ^^ 1 *!!^? 1 ^ 1 < C where c := c(k) is fixed. Then, there exists 



C fe71 



<PY t (£n+l -X n+ i). 



(Ixi.nl^ 1 + l^nl*- 1 )^ 2 

Remark 4.8 The condition in the Lemma means that the deviation £ n +i — x n+i has ex- 
actly the same order as the transport term t{\x\ )n \ k + |£i,n| fc ), up to a neglectable fluctuation 
corresponding to the variance of the Gaussian contribution in Yt . 

Proof. We assume w.l.o.g. that ^ n+ i - x n+ ± - ct(\x^ n \ k + |6,n| fe ) > and |£i >n | > |xi, n |. 
From (|4.19p we recall: 

PY t {in+l ~ X n+ i) = E[H t I Yt >Z n+1 -x n+1 ], 

H t := H} + H?:=^(D 7Yt ,DY t ) L2m + r Yt LY t . 

Recalling the chaos decomposition of LYt introduced in Proposition 14.21 see equation (|4.25j) . 
we get: 

PY^n+l ~ Xn+l) > E[H?I Yt >s n+1 - Xn+1 } ~ E[\H}\] > n hi9 " i - ,t) I Yt > (n+1 . Xn+1 ] 



lYt 



E 



lY t 



+ 



Ci 



m 



l,n 



using the bound for IE [ | ^f^/ 1 ] given by equation (|4.23j) . with C\ := C±(n,k, 1, if), in the last 
inequality. From equation (g2SD, there exists C 2 ■= C 2 (n,k), E[|X^ =2 ij(s>j_i(.,i))| 2 ] 1/2 < 



C 2 t^i,n\ k - 1 Et2(w^y * ^ t^ 1 ^ 2 ^- 1 , recalling \^ n \ > Kt 1 / 2 for the last 



< 
ft- 



inequality. Also h(go(.,t)) = £E[M \(u,t)]dW u + R where IE[| J^* I 2 ] 1 / 2 < ^l 2 \^ n \ k ~\ 
From (14.51). we write: 



Y t 



t i-t 

du\m(u, t, xi, n , £i,n)| fc + k du\m(u, t, xi tn ,^ n )\ k ~ 2 (m(u, t, xi >n , £i, n ), W" 
o Jo 

+Rt(x hn ,Zi,n) =■ (rn k + G k + i2?)(xi,n,6,n) =: m k + G k + R k , 



0,t\ 



for simplicity. Proposition 14.21 then yields: 



PY t (Cn+l -X n +l) > E 



fiE^u, t)]dW u 
lYt 

\ c 2 c 2 kt^ 2 \^ n \ k -^ 






\ K|^i, ft |2(fc-l)t3 ' ift 3 / 2 |6,n| fc - ■' 

Py tl i(£n+i -a^n+i) -ri(t,x,^). 
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From the martingale representation theorem and the above computations we identify G k 
/ *E[Mi(w,i)]dW u . Still from Proposition HJ we get: 



PYt (£n+l - x n+l) > E 



G*\ 



7r t 



J, 



G t fe +^>«n+i-^+i-mri^|<|G t fc |/2 
C 3 P[|^| > |G t fe |/2]V2 



+ ri(i,x,C) 



.(l^nl* -1 + l^nl*- 1 )* 3 / 2 
= Py t ,2(Cn+l -Xn+l) -r 2 (t,X,C), 

where G3 := C%{n, k). One easily gets that there exists c := c(k) > 0, m^ := m k (xi >n , £i,n) > 
ct(|xi in | fc + |£i, n | fc ).Thus, setting U^(x,C) ■= Cn+i - x n+ i - ct(\x ltn \ k + |6,n| fc ) and recalling 
as well that uf(x,£) > 0, one obtains that on the event {G k + R$ > U?(x,i), \R%\ < \G k \/2}, 
G k > 0. Hence: 

PYtHn+1 ~ X n +i) > H—lGt-^teUffaO^RZ^G^p] ~ r 2(t,X,^) 

/~<k 

> E[— lGfc> 2 [/ t fe (a;,O>0 I |i^|<Gj ; /2] ~~ r 2(^ x i 

> Ff G " II H U 1 C 3 F[ lYt >3M t ]^ 

[ 3M t G $> 2U tHx,O>0 l \R%\<G*/2 hy t <3Mt\ |^ lri |fe-l t 3/2 

-r 2 (i,x,£) 
- E [^T II G t fc >2Ct3/2(| a;iin |fe-i + |a,n| fc - 1 ) 1I |H t fc |<Gfc/2 II 7n< 3 - A/! «] - r 3(*,a;,0 

G-^GPfG^ > 2Gt 3 / 2 (|3;i i „| fc - 1 + |£i,n| fc-1 ), l^ fc | < G t fc /2] 



> 



3t 3 / 2 (|x lin |^i + l^nl*" 1 ) 

G 4 P[ 7 y t > 3M] 1 / 2 



t 3/2 |6,n 



fc-1 



?"3 (*,£,£)> 



(4.40) 



where we used that U k (x,£) < G£ 3 / 2 (|xi jn | fc_1 + |£i, n | fc_1 ) for the last but one inequality 
(compact sets of the metric) . The constant C is the one appearing in (|4,1U|) . To conclude it 
suffices to prove that 

P := nG^>2Ct^ 2 (\x 1 , n \ k - 1 + \^ n \ k - 1 ),\R^\<\G^\/2]}>C, (4.41) 

G 4 P[7y t > 3M t ] 1/2 , ^ C^CC /tA , 

M*,z,0 == AL lfc-1 +^3*,x,0< ^37^7] ufcZi— 77 life-n - ( 4 - 42 ) 



Indeed, plugging (|4.4ip and (|4.42p into (|4.40p gives the statement. Let us first prove (|4.4ip . 
Write: 

P > F[G k > 2C£ 3 / 2 (|x 1 , n | fc - 1 + ICi.nl^ 1 )] 

-P[G t fc > 2Gt 3 / 2 (K n | fc ~ 1 + iei.nl*- 1 ), \Rt\ > \G k \/2] 
> ¥[Af(0, 1) > 2G] - ¥[\R k \ > Ct*l 2 (\xi >n \ k - 1 + I6.nl* -1 )], G := C(n, k). 
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Thus, similarly to the proof of (|4.12p in Lemma 14.21 we can show that there exists C5 := 
C 5 {n,k) > 1 s.t. P[|i? 4 fc | > C^dx^nl*- 1 + |Ci,n| fc_1 )] <C 5 exp(-C^ 1 |a:i -" |2 + la " 12 ). Under 
the current assumptions, using standard controls on the Gaussian distribution function, this 
gives (|4.4ip for C := C(n, k) for K large enough. 

Recall now that \r A (t,x,£)\ < ,3/2,^-1 { Ul+c K ^ 2k + (C 3 + C 4 )F[ 7Yt > 3M t ]^ 2 + 
C3P[|i?f| > IG^I/2] 1 / 2 ) := X^i=i r 4j(^) x )^)- Under the current assumptions, we derive that 
for K large enough, r 4 i(£,x,£) < m z, 2{][X1 fprn^ nl*- 1 ) ' 0n the other hand ' writin g F 0- R * : | > 
|Gf|/2]V2 < Qp^fci > C ( | Xln |fe-i +|ei,„|*- 1 )t 3 /2]+P[|Gf| < C'dan,,,!*- 1 + l^i^l*- 1 )* 3 / 2 ]) 1 / 2 
we derive similarly to (|4.12p (see also the proof of the lower bound in Lemma I4.4J) that 
r4s(t, x,£) < -737271 — |fc-if|f — ffc=TT taking (7 small enough. Eventually, the same control 
holds true for r42(£,x,£), still from arguments similar to those used to derive (|4.12p . This 
concludes the proof. □ 

5 Potential Theory and PDEs 

In this section we are interested in proving Harnack inequalities for non-negative solutions to 

&u(z) = 0, z = (x, t) 6 M JV+1 , (5.1) 

with ££ defined in ()1.2p . Specifically, we consider any open set O C IR +1 , and any z £ 0, 
and we aim to show that there exists a compact K, <Z O and a positive constant Cjq such that 

supu < C/c u(z), (5-2) 

AC 

for every positive solution u to Jzfu = 0. We say that a set {^o, z%, . ■ . , Zf-} C is a Harnack 
chain of lenght k if 

u(zj) < Cj u(zj-i), for j = 1, . . . , fc, 

for every positive solution u of ££u = 0, so that we get 

u(z k )<C 1 C 2 ...C k u(z ). (5.3) 

In order to construct Harnack chains, and to have an explicit lower bound for the densities 
considered in this article, we will prove invariant Harnack inequalities w.r.t. a suitable Lie 
group structure. By exploiting the properties of homogeneity and translation invariance of 
the Lie group, we will find Harnack chains with the property that every Cj in (|5.3p agrees 
with the constant C/c in ()5.2p . As a consequence we find u{z k ) < C^u(zq), and the bound 
will depend only on the lenght of the Harnack chain connecting zq to z k . 

Let us now recall some basic notations concerning homogeneous Lie groups (we refer to 
the monograph [7] by Bonfiglioli, Lanconelli and Uguzzoni for an exhaustive treatment). Let 
o be a given group law on M^" 1 " 1 and suppose that the map (z, £) 1— > C _1 ° z is smooth. Then 
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G = (R , o) is called a Lie group. Moreover, G is said homogeneous if there exists a family 
of dilations {S\) x>0 which defines an automorphism of the group, i.e., 

S x ( z o C) = (5 x z) o (5x() . for a11 z > C e K^"*" 1 and A > 0. 
We also make the following assumption. 

[L] Jzf is Lie-invariant with respect to the Lie group G = (M iV+1 ,o, ((5a)a>o)> *- e - 

ij li, . . . , Y n and Z are left-invariant with respect to the composition law of G, i.e. 

Y j (u((o-)) = (Y j u)((o-), j = l,...,n, 
Z(u((o.)) = (Zu)((o.), 

for every function u G C 00 ^ 4 " 1 ), and for any ( G M Af+1 ; 

wj Y"i , . . . , Y n are (^-homogeneous of degree one and Z is (^-homogeneous of degree 
two: 

Yj (u (S\z)) = A (Yju) (5\z) , j = 1, . . . , n, 
Z(u(5 x z)) = X 2 (Zu)(5 x z), 

for every function u G C°°(R Ar+1 ), and for any z G M iV+1 , A > 0. 

To illustrate Property [L] we recall the Lie group structure of the Kolmogorov operator 
corresponding to k = 1 in (jl.4|) . 

Example 5.1 (KOLMOGOROV OPERATORS) Jzf := ^A Xln + Ya=X x $xi ~ ®t- Tfie kol- 
mogorov group is K = (M n+2 ,o, <J A ) ; where 

n 

(x, t) o (£, T ) = (xi >n + 6,n, £n+i + Cn+i - ^ »iT, t + t) , 5 A (x, i) = (Axi, n , A 3 x n +i, A 2 t) . 

Clearly, Jzf can 6e written as in (|1.2p lozf/i Yi = d Xi , i G [l,rt], and Z = X]iLi Xid Xn+1 — dt, 
and satisfies [L] . 

It is known that the composition law o is always a sum with respect to the t variable (see 
Propostion 10.2 in [23]). Moreover, the family (#a)a>o ac ^ s on ^ Ar+1 as follows: 



8\ (x\,X2, ■ ■ ■ ,XN,t) = (X ai xi,X (72 X2, ■ ■ ■ , \ aN X]\T, \ 2 t) , for every (x,t) G 



pAf+l 



where a = (<Jx,o~2, ■ ■ ■ , o"jv) £ f^ is a multi-index. The natural number Q = X^fe=i °~k + 2 
is called the homogeneous dimension of G with respect to 5\. We shall assume that Q > 3. 
Observe that the diagonal decay of the heat kernel on the homogeneous Lie group is given 
by the characteristic time scale t~^~ 2 '' 2 . For the above example we have Q = n + 3 + 2, 
matching the diagonal exponent in (|1.5p (Q — 2)/2 = (n + 3)/2. 
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Write the operator «£? as follows 

N N 

% = ^ Oijitfdxuxj + / J b j (x)d Xi - d t , 

for suitable smooth coefficients a^j's and bj's only depending on the vector fields Yq, ... ,Y n . 
As n < N, «£? is strictly degenerate, since the rank(j4(x)) < n at every x (here A(x) := 
(aij{x))- ■ w-, ,). In Example 15.11 we see that rank(A) never vanishes. We say that J2? is not 
totally degenerate if 

[B] for every xeR^ there exists v G R N \ {0} such that (A(x)u, v) > 0. 

This property holds for a more general class of operators. Indeed, if J§f satisfies [H] and [L], 
then there exists a v G R \ {0} such that 

(A(z)i/, z/) > 0, for every x G 1^. (5.4) 

We refer to Section 1.3 in the monograph [7j for the proof of this statement. 

Fix now T > and define / := [0,T]. We call diffusion trajectory any absolutely contin- 
uous curve on / such that 

11 
l'{s) = ^2uj k (s)Y k (j(s)), for every s G I, (5.5) 

fe=l 

where cji, . . . ,u n are piecewise constant real functions. A drift trajectory is any positively 
oriented integral curve of Z. We say that a curve 7 : [0, T] — > M Af+1 is J£ -admissible if it is 
absolutely continuous and is a sum of a finite number of diffusion and drift trajectories. 

Let O be any open subset of IR^" 1 " 1 , and let zq G O. We define the attainable setsrf Zo := A Zo 
as the closure in O of the following set 

A Zo = {z G O : there exists an Jz? -admissible path 

7 : [0,T] -> such that 7(0) = z ,7(T) = z}. 

The main result of the section is the following 

Theorem 5.2 Let S£ be an operator in the form (|5.ip satisfying [H] and [L], Ze£ C R^" 1 " 1 
be an open set, and let zq G O. Then, 

for every compact set K, C Int (<e4 ) , supn < Cx u(^o), (5-7) 

K 

for any non-negative solutions u to ^£u = in O. Here C/c is a positive constant depending 
on O, /C, zq and on ££ . 

We recall that a Harnack inequality for operators satisfying [H] and [B] is due to Bony 
(see [9]). Another result analogous to Theorem 15.21 is given in |15|. Theorem 1.1] by Cinti, 
Nystrom and Polidoro, assuming [L] and the following controllability condition: 
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[C] for every (x,t), (£, r) G M iV+1 with t > t, there exists an Jzf-admissible path 7 : [0, T] 
m such that 7(0) = (x,t), j(T) = (f,r). 



Our Theorem 15.21 improves Bony's one in that it gives an explicit geometric description of 
the set K, in (|5.7|) . Also, it is more general than the one in [T5], since [L] and [C] imply [H] 
(see Proposition 10.1 in |24|). 

The proof of Theorem 15.21 is based on a general result from Potential Theory. In Section 
15.11 we recall the basic results of Potential Theory needed in our work, then we apply them 
to operators Jz? satifying [H] and [L]. We explicitly remark that condition [L] is not satisfied 
by the Kolmogorov operators (|2.7|) and (|2.8j) . 

5.1 Potential Theory 

For the first part of the section, we assume Jz? to be a general abstract parabolic differential 
operator satisfying [B] and [L]. 

Let O be any open subset of M Ar+1 . If u : O — V R is a smoothfunction such that Jzfti = 
in O, we say that u is ^-harmonic in O. We denote by %{0) the linear space of functions 
which are Jz? -harmonic in O. 

Let V be a bounded open subset of R Ar+1 with Lipschitz-continuous boundary. We say 
that V is Jzf '-regular if, for every zq G dV, there exists a neighborhood U of zq and a smooth 
function w : U —> M. satisfying 

io(>o) = 0, Sfw{z )<0, w>0inFn[/\{z }. 

Note that the function ip(x,t) = i + - arctant verifies 

< ip < 1, ^V < in l^ 1 . (5.8) 



As a first consequence of (|5.8p . the classical Picone's maximum principle holds on any bounded 
open set O C M Ar+1 . Precisely, if u € C 2 {0) satisfies 

S£u > in O, limsupu(2;) < for every (" G dO, 

then it < in O (see e.g. Bonfiglioli and Uguzzoni [8j). Then, for every Jzf-regular open set 
V C IR^" 1-1 , and for any ip G C(5V) there exists a unique function H^ satisfying 

HX EH(V), lim HX(z) = ip(() for every C G SV. (5.9) 

Moreover, IlY, > whenever 99 > (see Bauer [1] and Constantinescu and Cornea |18|). 
Hence, if V is Jzf-regular, for every fixed z G V the map 99 1— > H^(z) defines a linear positive 
functional on C(dV,M). Thus, the Riesz representation theorem implies that there exists a 
Radon measure /u]f, supported in dV, such that 

HY(z)= I (piQdfiYiQ, for every <peC(dV,R). (5.10) 

Jav 
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We will refer to \x], as the ££ -harmonic measure defined with respect to V and z. 

A lower semi-continuous function u : O — >] — oo, oo] is said to be Jif-superharmonic in O 
if u < oo in a dense subset of O and if 

u(z)> [ u(()d^(C), 
JdV 

for every open Jzf-regular set V C V C O and for every z E V. We denote by S{0) the 
set of Jjf-superharmonic functions in O, and by S (O) the set of the functions in S{0) 
which are non-negative. A function v : O — > [—00, oof is said to be J£-subharmonic in O if 
—v £ S(0) and we write S(0) := —S(0). Since the collection of £-regular sets is a basis for 
the Euclidean topology (as we will see in a moment), we have S{0) n«S(C) = H(0). 

This last property and Picone's maximum principle are the main tools in order to show 
the following criterion of Jzf-superharmonicity for functions of class C 2 (a proof can be found 
in the monograph [TJ Proposition 7.2.5]). 

Remark 5.3 Let u E C 2 (0). Then u is J£ -superharmonic if and only if J£u < in O. 

With the terminology of Potential Theory (we refer to the monographs [H [18] ) , the map 
jj-W+l 3 q ,_). "^(C) is said harmonic sheaf and (JL N+1 ,1-L) is said harmonic space. Since 
the constant functions are Jzf -harmonic, the last statement is a consequence of the following 
properties: 

- the Jz? -regular sets form a basis for the Euclidean topology (by (|5,4p , ££ is a not totally 
degenerate operator, so that this statement is a consequence of [9[ Corollaire 5.2]); 

- % satisfies the Doob convergence property, i.e., the pointwise limit u of any increasing 
sequence {u n } n of Jzf-harmonic functions, on any open set V, is Jz?-harmonic whenever 
u is finite in a dense set T C V (as in [241 Proposition 7.4], we can rely on the weak 
Harnack inequality due to Bony stated in Theoreme 7.1]); 

- the family 5(M +1 ) separates the points of R +1 , i.e., for every z, ( € R +1 , z / (, 
there exists u € 5(R JV+1 ) such that u(z) 7^ u(C)- 

This last separation property is proved in Lemma 15.51 We will in fact show a stronger 
result: actually, the family S (R w+1 ) n C(R Ar+1 ) separates the points of R^"*" 1 . A harmonic 
space (R +1 ,7i) satisfying this property is said to be a ^-harmonic space. 

In order to prove the separation property we use a fundamental solution T of Jz? . To 
prove the existence of a fundamental solution we now rely on condition [H] that we assumed 
to be in force through the paper. We recall that a fundamental solution is a function F with 
the following properties: 

i) the map (z,C) ^ T(z,£) is defined, non-negative and smooth away from the set {(z,Q G 



oJV+l 



x 



l 'z*Q\ 



a 



',) for any z G R^" 1 " 1 , T(-, z) and T(z, •) are locally integrable; 
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:) for every <j> G C{f(R N+1 ) and z G R N+l we have 



/RJV+l 



00(C) dC = / r(z, C W(C) <*C = -#*); 



iv) Jz?r(-,£) = — <5^ (Dirac measure supported at (); 

v) if we define T*(z, £) := r(£, z), then T* is the fundamental solution for the formal adjoint 
Jzf* of ££ , satisfying the dual statements of Hi), iv); 

vi) r(x,t,£,r) =0 if t<r. 

Remark 5.4 Assumption [H] implies the existence of a smooth density p{t,£,x)dx := P^f-X* G 
dx], i > 0, for the process (Xt)t>o associated to S£ see e.g. Stroock |37^ or Nualart f3^ . Ac- 
tually, 

r(M,£,r) :=p(t-r,^,x) 

is a fundamental solution for ^f in the above sense. Indeed p satisfies the Kolmogorov equation 
Jzfp = 0, in R N+1 \{(^, t)}. We re/er to Bonfiglioli and Lanconelli JBj for a purely analytic 
proof of existence of fundamental solutions for operators satisfying [H] , [L] . 

If condition [L] holds, then we also have: 
vii) T(z, C) = r(a oz,ao() for every a, z, C € M^" 1 " 1 , z ^ (; 

viii) T(S x (z), <5 A (0) = A-« +2 r(z, c), «, C e ^ N+1 , z± C, A > o. 

We next prove the separation property for Jz? by adapting the argument in [14\ Proposition 
7-1]. 

Lemma 5.5 For every Z\,z% G K^" 1 " 1 , Z\ ^ Z2, there exists a function u G S (R N+l ) n 
I +1 ) sitc/i £/ta£ tt(zi) 7^ u(/22). 



Proof. Let us denote ^ = (xi,U) for i = 1,2. First we suppose that ii < ti- The properties 
of T yield that there exists zq = (xo,to) with to > such that T(zq,0) > 0. On the other 
hand, since [H] and [L] yield [B], there exists a «£?-regular open set Vo containing the origin, 
a small r$ > and a large Ao > 1 such that 

U ro CV C6 Xo {U ro ), U ro ={{x 1 ,...,x N ,t)eR N+1 :\x i \ < r ,\t\ <r }. (5.11) 

By the smoothness of T, there exists e > such that T > in the set zq o U e . For a fixed 



A G 

set 



0, 



fa-fa 



' V 2 (*0+e) 



and a non-negative function 93 G C^°(z2 o {8\(zq o U £ )) n {i < £2}), we 



u v (z)= T(z,()p(()d(, zeR N+L . (5.12) 



Hence, we obtain u^ G C 00 ^^" 1 " 1 ), u 9 > and Jzfw^ = —99 < 0, so that, by Remark [53 
-u^ G S(M. N+1 ). Moreover the choice of ip implies that u v {z\) = and 14^(^2) > 0. 
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In the case t± = t2, X\ 7^ X2, we consider the sequence 

O n (z 2 ) = {(£R N+1 :r(z 2 ,()>n Q - 2 }, n G N. (5.13) 

We note that O n {z2) shrinks to {z 2 } as n — > 00, by property viii) of the fundamental solution. 
For any cp n G C™(O n (z2)) such that J cp n = 1 and (p n > 0, we define u Vn as in (j5.12|) . Then, 
Uip n is a smooth non-negative function in IR^" 1 " 1 satisfying 5£u [pn < 0, and so u Vn is Jzf- 
superharmonic. It holds 

r(z 2 ,c)^n(CMC>7i Q - 2 

hi 



u Vn (z 2 ) = / r(z 2 ,C) <Pn(Q dQ>n Q for every n G N; 

••+1 

%> n (^l) < max r(zi,C) = C, 



C6Cl(^2) 

where C is a real positive constant independent of n. This ends the proof. □ 

We summarize the above facts in the following 

Proposition 5.6 Let J£ be an operator in the form ()5.ip and assume that [H] and [L] are 

satisfied. The map Ti which associates any open set O C M-" +1 with the linear space of the 
^-harmonic functions in O is a harmonic sheaf, and (R + , T-L) is a QS-harmonic space. 

A remarkable feature of a QS-harmonic space is that the Wiener resolutivity theorem holds 
(see [Hdl]). In order to state it, we introduce some additional notations. We recall that if 
O C M. +1 is a bounded open set, then an extended real function / : dO — > [—00, 00] is called 
resolutive if 

infZZy = supWf =: Hf G rl(O), 



where 



U° f = {u 6 S{0) : infu > -00 and liminf u(z) > /(C), VC € dO}, 



Wf = {u€ S{0) : supn < 00 and limsupu(z) < /(C), VC G dO}. 

We say that Hf is the generalized solution in the sense of Perron- Wiener- Brelot to the 
problem 

ueH(O), u = f on dO. 

The Wiener resolutivity theorem yields that any / G C(dO,M) is resolutive. The map 
C(dO,M.) 9/i—)- Hf(z) defines a linear positive functional for every z G O. Again, there 
exists a Radon measure /if on dO such that 

Hf(z)= [ /(C)d^(C). (5.14) 

ho 

We call /if the J£ -harmonic measure relative to O and z, and when O is ££ -regular this 
definition coincides with the one in (|5.10p . Finally, a point C G dO is called J£ -regular for O 
if 

lim Hf(z) = f(Q, for every feC(dO,R). (5.15) 

Obviously, O is .if-regular if and only if every C G dO is Jzf-regular. 
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5.2 Harnack inequalities 

Let O C IR^" 1 " 1 be an open set. A closed subset F of O is called an absorbent set if, for any 
z G F and any Jzf -regular neighborhood V C V C of z, it holds //]f (SV^ \ F) = 0. For any 
given zq G O we set 

j^so = {F C O : F B zq, F is an absorbent set}. 

Then, 

O zo = f| F (5.16) 

is the smallest absorbent set containing zo. The Potential Theory provides us with the 
following Harnack inequality. Let (R N+1 ,H) be a ^-harmonic space, let O be an open subset 
ofR N+1 and let z G O. Then, 

for every compact set K C Int (O zo ) , supu < Ck u(zq), (5-17) 

K 

for any non-negative function u G H((D). Here Ck is a positive constant depending on 
0,K,zq. We refer to Theorem 1.4.4 in [3] and Proposition 6.1.5 in [18]. Proposition 15.61 
implies that (|5.17p applies to our operator ^f. We summarize the above argument in the 
following 

Proposition 5.7 Let Jzf be an operator in the form (|5.ip satisfying [H] and \L], let O C 
K +1 be an open set, and let zq G O. Then, 

for every compact set K C Int (0 ZQ ) , supu < Ck u{zq), 

k 

for any non-negative solutions u to Jzfu = in O. Here Ck is a positive constant depending 
on O, K, zq and on jSf. 

In order to prove Theorem 15.21 we give the following 

Lemma 5.8 Let Jzf be an operator as in (|5.ip satisfying [H] and [L], and let O be an open 
subset o/M +1 . For any given zq £ O, we have s^ ZQ C O zo with £/ zo defined in (|5.6[) . 

Proof. Since O zo is a closed set, and &/ Zo is the closure of the set A ZQ defined in (|5.6p . it is 
sufficient to show that A zo C O zo . By contradiction, assume that ~z G ^2; \ O zo . Then, there 
exists an Jzf-admissible path 7 : [0, T] — > O such that 7(0) = zq, 7(F) = z. 
We set 

t i: =inf{t>O: 7 (]t,F])nO 2o =0}. 

Note that, since O \ O zo is an open set containing z and 7 is a continuous curve, there 
exists an open neighborhood U C O of z such that (7 n O zo = 0, and a positive a satisfying 
7(]F— a, T\) C f7. Hence, £1 G [0, F[ is well defined and we have 7(4) ^ O zo for every t G]£i, F]. 
Again, by the continuity of 7, we have 

zi =7(*i) eO Z0 . 
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Let V C V C O be a J2?-regular neighborhood of z\ with z ^ V. Arguing as above, we can 
find £2 ^li^I such that 7([£i,i2[) C V and z<i = 7(^2) £ 9V. Consider any neighborhood W 
of Z2, such that WcO\ O zo . Let <p £ C(dV) be any non-negative function, supported in 
WndV, and such that 9? (22) > 0. Recalling that the harmonic function H^ is non-negative, 
we aim to show that 

H^izx) > 0. (5.18) 

By contradiction, we suppose that nX vanishes at z\. In other terms, HX attains its minimum 

value at z%, then Bony's minimum principle implies HX = in 7([ii,i2[)- As a consequence, 

since HX satisfies (|5.9p . 

hmHY( 7 (t)) = Q. (5.19) 

t-rt 2 

On the other hand, by the choice of <*p 



lim F;(z) = V 9(z 2 )>0. 



This contradicts (|5.19|) and proves (|5.18|1 . By using representation (|5.10|1 of HX. in terms of 
the j£f -harmonic measure, (|5.18p reads as follows 

HY(zi)= [ ¥>(C)<K(C) >0, then ^(dV DW) > 0. (5.20) 



On the other hand, z\ belongs to the absorbent set O zo , so that fiY (dV \ O zo ) = 0. But this 
clashes with (|5.20p . being W O O \ O zo . This accomplishes the proof. □ 

Proof of Theorem \5.2l It is a plain consequence of Proposition 15.71 and Lemma 15.81 □ 

As the following proposition shows, we are able to give a complete characterization of the 
set O zo if g/ Zo is an absorbent set as well. 

Proposition 5.9 Let «$? be an operator as in (|5.ip satisfying [H] and [L], let O C M" +1 6e 
an open set, and let zq £ O. If £/ ZQ is an absorbent set, then £/ ZQ = O zo . 

Proof. The claim directly follows from Lemma 15.81 recalling the definition of O zo . □ 

The first statement in next proposition is a classical result in abstract potential theory 
(see e.g. [H Theorem 1.4.1] and |18t Proposition 6.1.1]). For the convenience of the reader, 
we explicitly give here its simple proof. 

Proposition 5.10 Let 5£ be an operator as in (|5.ip satisfying [H] and [L], let O C M" +1 
be an open set, and let Zq € O. Assume that there exists a solution u > to J£u = in O 
such that u = in £/ ZQ and u > in 0\ g/ zo . Then g/ Z() is an absorbent set, and g/ zo = O zo . 

Proof. Since u is continuous and non-negative, 

£f zo ={zeO: u{z) < 0} 

is a closed subset of O. Let z £ £/ zo , and let V C V C O be a S£ -regular neighborhood of z. 
As u e W(0), we have 

> u (z) = / u(C) d^(() > 0, so that /J,X(dV \^ Zo ) = 0. 

JdV 

Hence £/ Zo is an absorbent set. The last statement plainly follows from Proposition 15. 9[ □ 
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5.3 Lifting and Harnack inequalities 

We first consider the PDE (|2.7|) for k = 2. Note that, in this case, it is equivalent to 
and reads as follows 

^ = l^ n + \xi,n\ 2 d Xn+1 -d t . (5.21) 

It is homogeneous with respect to the following dilation 

6 X (x, t) = (Axi, n , A 4 x n+ i, X 2 t) . (5.22) 

Even if Jz? does not satisfy [L]-i), it has a fundamental solution T which shares several 
properties of the usual heat kernels. We remark that, since Jz? does not satisfy the control- 
lability condition [C], the support of F is strictly contained in the half space {t < r}. 

We next show that Jz? can be lifted to a suitable operator Jz? in the form (|5.ip satisfying 
both [H] and [L]. By adding a new variable y = y\ iTl G M n , we define the following vector 
fields on R 2n+2 

n 

Yi = Yi = d Xi ,i£ll,nl Z = \x 1 , n \ 2 d Xn+1 +Y / x i d yi -d t . (5.23) 

Clearly, if we denote v(x,y,t) = u(x,t) for any u S C 00 (M n+2 ), we have 

Yiv(x,y,t) = Yiu(x,t), Vie[l,nJ, Zv(x,y,t) = Zu(x,t), 

then, if we consider the lifted operator Jz? = | Y17=l ^i + %■> we ^ no - -^v{x, U, t) = ^fu{x, t). 
By a standard procedure (see e.g., Chapter 1]), we explicitly write the group law o of the 
homogeneous Lie group G = (lR 2n+2 , o, (5\)\ > o) such that Jz? is G-Lie-invariant: 

(x,y,t)o(£,r],T) = (xi,n+£l,n, #n+l +£n+l + 2(xi, n , 171^) - t\x 1jH \ 2 , y 1>n + T] 1 ,n ~ TXl >n , t + t), 

(5.24) 
and the dilation 6\: 

5 x (x, y, t) = (Xx 1>n , A 4 a; n+ i, A 3 yi,n, \ 2 t) . (5.25) 

Therefore, the lifted operator Jz? satisfies [H] and [L].In the sequel we will consider admissible 
paths in the following form 

n 

j'(s) = 5^^(7(8)) + Z(j(s)), s € [0,r], 

3=1 

for some constant vector u = (u)\, . . . ,u n ),^/(0) = (x,y,t). Its explicit expression is 

/ s 3 s 2 \ 

j(s) = f xi,„ + sw, x n+ i + s|xi in | 2 + s 2 (xi jn ,w) + — |w| 2 ,y + sxi, n + — w,i-sj . (5.26) 

In order to prove an invariant Harnack inequality for the non-negative solutions to Jzfv = 
0, we describe the sets O zo and g/ zo in the case when zq is the origin and 



O 



{{x,y,t)eR 2n+2 \ \x lyn \ <l,-Kx n+l <l,\y\ <l,-l<t<l}. (5.27) 
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Lemma 5.11 Let O be the open set defined in (|5.27p . and let zq = (0,0,0). Then 

£/ Z() = {(x,y,t) e O | < x n+ i < -t,\y\ 2 < -tx n+ i} , (5.1 

and O, 



'ZO 



Proof. In order to prove (|5.28p . we consider any Jzf- admissible curve 7 in O. In our set- 
ting, the components Xn+ijJ/i,™ and t of every diffusion trajectory are constant functions. 
Moreover, any drift trajectory 7 : [0, T] — >• O starting from (x,y,T) is given by 

7(5) = (x 1 , n ,x n+1 + s\x 1<n \ 2 ,y + sxi tTl ,t-s). (5.29) 

Hence, any Jzf-admissible curve 7 : [0,T] — > O with 7(0) = (0,0,0) is given by 



7(a) = xi 



M, / ^2\c k \ 2 h k (r)dr, / ^2c k I Ik (r)dr,- ^h k (r)dr , sG[0,T]. 
- 70 fc =i - 70 fc=i - 70 fc=i / 



Here I\,...,I m are disjoint intervals contained in [0, T] and Ij h denotes the characteristic 
function of I k . The function x\ iH is constant on every I k , and any c k is a constant vector 
such that |cfc| < 1 for A; = 1, . . . , m. As a consequence of the Holder inequality we find 

<s4 ^ {(^yj*) € O I < x n+ i < -t, |y| 2 < -tx n+ i} . 

In order to prove the opposite inclusion, we consider any point 

(x,y,t) £ {(x,y,t) £ O \ < x n+ i < -t,\y\ 2 < -tx n+1 ) , y ^ 0, 

and we show that there exists a Jzf- admissible curve 7 = 71 + 72 + • • • + 75 contained in O, 
which steers (0,0,0) to (x,y,t). To this aim, we fix a small positive e, that will be specified 
in the sequel, and we set 

-tx n+ i - \y\ 2 



£ x n+1 -2\y\(l-e)-t(l-e) 2 ' 

Note that -tx n+1 + 2\y\t(l - e) +t(l - e) 2 > (\y\ + t(l - e)) 2 , so that < s £ < -t. We set 
Xl >n = ^pry and we choose 71 as a diffusion trajectory connecting (0,0,0) to (5?i jn , 0,0,0), 

and 72 : [0, s £ ] — > M. 2n+2 as a drift trajectory starting from (xi in , 0,0,0). Hence, according 
to (|5.29|) . we find 72(s £ ) = (x\^ n ,s £ {l — e) 2 ,s £ ^^y,—s £ ). Then, by a diffusion trajectory 

73, we connect 72(s £ ) to the point ( 7^- \\-\ Vi s e(l — e ) 2 ; s ei=fy) ~ s e ) • We next consider a 
drift path 74 : [0, —t — s £ ] — > M. 2n+2 which, by (|5.29p . and by our choice of s E , steers the end 
point of 73 to ( (_t_ )\-\ y^n+i,y,t) ■ Finally, we can find a diffusion path 75 connecting 
7 4 (-t-s £ ) to (x,y,t). _ 

Clearly, 7 = 71 + 72 + • • • + 75 is a Jzf-admissible curve of R 2n+2 connecting (0, 0, 0) to 
(x,y,t). Next we prove that, for sufficiently small e, the trajectory 7 is contained in O. To 
this aim, as the set O is convex and the paths 71,72, • • • ,75 are segments, we only need to 
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show that the end-points of 71,72,73,74 belong to O. The inequalities — 1 < ' y ' J-_ — — < 1 
directly follow from the definition of s £ , for sufficiently small positive s. The other inequalities 
are a plain consequence of the fact that < s £ < — t < 1, as previously noticed. Since g/ ZQ is 
the closure of the set of the points that can be reached by a Jzf-admissible path, we get 



{(x,y,t) E O I < x n+ i < -t,\y\ 2 < -tx n+1 ] C 



*z - 



This concludes the proof of (|5,28p . 

To complete the proof, by Proposition 15.101 it is sufficient to find a non-negative solution 
v of S£v = 0, such that v = in £/ ZQ , and v > in O \ £/ ZQ . Let ip be any function in C(dO), 
such that ip = in dO n g/ zo and ip > in <90 \ ^ . Then the Perron- Wiener-Brelot solution 



v := i?~; of the following Cauchy-Dirichlet problem 



I ££ v = in O 

\v = ip in 90 

is non-negative. Next we prove that v > in 0\s^ ZQ . By contradiction, let (x,y,t) E O\g/ Z0 
be such that v(x, y, i) = 0. Then (x, y, t) is a minimum for u, so that from Bony's minimum 
principle Theoreme 3.2] it follows that v(xi tn ,x n+ i,y, t) = <p(xi <n ,x n+ i,y,t) = 0, for 
every x\ :n £ d(] — 1, l[ n ). Since every point (xi :n ,x n+ i,y, t) is regular for the Dirichlet 
problem, and belongs to dO\.s/ Z0 , we find a contradiction with our assumption on ip. Suppose 
now that there exists (x,y,t) E g/ Zo such that v(x,y,t) > 0. Since every point of the set 
dO fl g/ Zo is J2?-regular, v is continuous in srf ZQ . Hence there exists a (x,y,t) E £/ Z() such 
that v(x,y,t) = max^ v > 0. By Bony's minimum principle we have v(xi jn ,x n+ i,y, t) = 
(p(xi jn ,x n+ i,y,t) > 0, for any x\^ n E d(] — 1, l[ n ), and this fact contradicts our assumption 
on 93. □ 

Next we introduce some notations to state a Harnack inequality which is invariant with 
respect to the group law o defined in (|5.24p and the dilation 5 r introduced in (|5.25p . Consider 
the box Q r =] — r,r[ n x] — r A ,r 4 [x] — r 3 ,r 3 [ n x] — r 2 ,0], and note that Q r = 5 r Q\. For every 
compact set K, C Qi, for any positive r and for any zq E M 2n+2 we denote by 

Qr(zo) = z o 8 r Qi = {z o 5 r ( I C G Qi}, fc r (za) = zqo 5 r X. (5.30) 

Corollary 5.12 For every compact set K, C \(x,y,t) E Q\ \ < x n+ \ < —t, \y\ 2 < — tx n+ \} , 

r > and z$ E M. 2n+2 there exists a positive constant C;c, depending only on Jz? and /C, such 
that 

sup v < C/cv(z ), 

for every non-negative solution v of Jifv = on any open set containing Q r (zo). 

Proof. Consider the function w(z) = v(zo o 5 r z) . By the invariance with respect to 5 r and o, 
we have ££w = in Q\. Aiming to apply Theorem 15.21 we consider the open set O defined 
in (|5.27p . and we note that O n {t < 0} C Q\. Then w is defined as a continuous function 
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on dO n {t < 0}. We extend wtoa continuous function on dO, and we solve the boundary 
value problem Jzfu; = in Qx, with w = w in dO. Then we apply Theorem 15.21 and Lemma 
15.111 and we get sup^ w < Cjcw(0,0,0). By the comparison principle we have w = w in 
On {t < 0}, then the claim plainly follows from the inclusion /CcOn{t<0}. □ 

We are now ready to build a Harnack chain for (|5.2ip by using the following set 



K 



{(x,y,t)em. 2n+2 \ |xi,n| <\^<x n+l <\,\ y \ <§,* = -£} (5-31) 



which is a compact subset of {(x,y,i) G Q\ | < x„,+i < —t, \y\ 2 < — ix n +i}. Before doing 
that for k = 2 only, we extend the above procedure to equations (|2.7f) and (|2.8p for k > 2. 

We next show that, in both cases ([2.7|) and (|2.8|) . j£f can be lifted to a suitable operator 
«Sf in the form (|5.ip satisfying [H] and [L]. We introduce a new variable y G R( fe-1 ) n , that will 
be denoted as follows y = (y x , y 2 , • • • , 2/(&-i)), with y,- = (y^i, . . . , y jn ) G R n for j G [1, k - 1]. 
We then define the lifted vector fields on R kn+2 : 

ft— 1 n 

?i = I* = 3*, * € [1, n], Z = Z + Y, S x i 5 ^ ; ^ 5 - 32 ) 

i=l i=l 

where Z = |xi, n | fc ^„ +1 - <9 4 for <|2Z}, and Z = £" =1 xjd x „ +1 " dt for d^D" If we denote 
v(x,y,t) = u(x,t) for any u G C°°(IR n+2 ), we have 

Yiv(x,y,t) = Yiu(x,t), ViG[l,ra], Zv(x,y,t) = Zu(x,t). 

Then, setting Jz? = ^ Ys?=l ^i + ^> we plainly find J£v{x, y, t) = Jt?u(x, t). 

Since dim(Lie{Y"i, . . . , y n , Z}) = /en + 2 and rank(Lie{Y"i, . . . , Y n , Z}(x, y, £)) = kn + 2 
at every point (x,y,t) G IR +2 , Theorem 1.1 in [6] yields the existence of a homogeneous 
Lie group G = (K fcn+2 ,o, (<5a)a>o) such that Jzf is Lie-invariant on G. Therefore, the lifted 
operators Jz? satisfy [H] and [L]. The dilation 5\ acts as follows: 

S x (x,y,t) = (\x 1:n ,\ k+2 x n+1 ,\ 3 y 1 ,...,\ k+1 y k ^ 1 ,\ 2 t^ , (5.33) 

for every (x,y,t) G M. kn+2 , and A > 0. We next aim to apply Theorem 15.21 in order to 
prove a Harnack inequality on the lifted space M. +2 . For any to G L 2 ([— T, T],M n ) for every 
0, y, t) G M fcn+2 and T > 0, we denote by 7 : [-T, T] -»■ M fcn+2 the solution of the Cauchy 
problem 

\l'{s) = Y:; =l ^{s)YMs)) + Z{l{s)), se[-T,T], 

\j(0) = (x,y,t). 

In order to simplify the notation, in the sequel we will denote the solution of (|5.34p as 

l(s) = (xi,n{s),x n+ i(s),y(s),t(s)), sE[-T,T]. (5.35) 

Note that t(s) = t — s for every s G [-T, T], so that t(T) = t — T. 
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The composition law "o" of G is related to (|5.34p as follows: if (x,y,i) = y(T) is the 
end point of the path 7 defined by (|5.34p with 7(0) = (0, 0, 0) and (x, y,t) = j(T) is the end 
point of the path 7 defined by (|5.34|) with 7(0) = (£, rj, r), then 

ix,y,t) = (£,tj,t)o (x,y,t) , (5.36) 

with t = —T (see for instance Corollary 1.2.24 in [7]). The above identity also holds when 
computing 7 at s = —T. In particular, if we choose any U G IR n , and t > 0, we let T = 
t,uj(s) = uJ for any s G [— t, t\, we find 

x= (-iu, -T^\u}\ k ^J , t = T + t, 

( /* - k \ (5 - 37) 

» = Kl,n ~ tu, Cn+l ~ / |6,n~sw| ds J . 

According with Remark I5.4|, the fundamental solution r of Jz? exists and is invariant with 
respect to the group operations (|5.36p and (|5.33p . Then function 

T(x,t,Z,T)= I r(x,y,t,£,0,T)dy (5.38) 

JR(fc-l)" 

is a fundamental solution to Jz? and gets from T the following invariance properties: 

T{x,t^,r)=T(x,t,0,0), 

1 ( 5 - 39 ) 

r(Axi jn ,A +2 x n+ i,A 2 t,0i jn+ i,0) = An+fc+2 T(x,t, 0i, w +i , 0) , 

for every (£,r)(a;,£) 6 R n+2 ,A > 0, where (x,t) is defined in ()5.36p . In the following remark 
we summarize the above properties when (x,t) has the form (|5.37p . 

Remark 5.13 For every (£,t) G R n+2 ,£ > and for any constant vector uj € M n , we have 

r Kl.n - VtUJ, Cn+l ~ / £l,n--7=U ds,T + t,£,,TJ = n+k — T^ - U, -j^, 0l,n+l, °) • 

We next focus on the attainable set ^4 of the unit cylinder 

= {(i,i,,{)eR k+2 |xi, n | <l,-l<x„+i<l,|y| < 1,-1 < i < l}, (5.40) 

with respect to the point zq = (0,0,0). Here |xi,n| an d \y\ denote, respectively, the Euclidean 
norm of the vectors x\ >n £ R n and y G M.( > n . 

Unlike the case k = 2, as k > 2 we are not able to give a complete characterization of the 
sets <s/ zo and O zo as we did in Lemma 15.111 We will consider instead the differential of the 
end point map related to (|5.34p to find some interior points of <e^ . With obvious meaning 
of the notations, we set (x(T), y{T), t(T)) = j(T), we note that t(T) = t — T, and we define 

E:L 2 ([0,T])^R kn+1 , E(u) = E(u,x,y,t,T) := (x(T),y(T)). (5.41) 
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We refer to the classical literature (see e.g. \i2\ Theorem 3.2.6]) for the differentiability 
properties of E. We next show that the differential DE(lo) of E, computed at some given 
u) € L 2 ([0,T]) is surjective. Hence E(u) is an interior point of ^4 , so that we can apply 
Theorem [ 



Lemma 5.14 Let w be any given vector ofW 1 such that uij 7^ for every j 6 [1, n] . Consider 
the solution 7 to the problem (|5.34|) . with uj = w. Then DE(uj) is surjective. 

Proof. By the invariance of the vector fields Y{, i £ [1, n], and Z with respect to the homoge- 
neous Lie group G, is not restrictive to assume (x,y,t) = (0,0,0) and T = 1. To prove our 
claim, we compute 

DE(uj)u = lim - (E(u + hu) - E(u)) , 
h—*0 h 

where 

u)(s) = v for s E [a, b], a,b £ [0, 1], a < b, v is any vector of W 1 , 

£(s) = for s<£[a,b}. (5.42) 

In the sequel, we denote by ^ h (s) = (x h (s) , y h (s) , t h (s)) the solution of (|5.34p relevant to 
w + hui. Clearly, t h {s) = — s, and x h (l) =w + hv, so that 

limi(4 in (l)-x lin (l))= V . (5.43) 

We next show that, for every j £ [1, n] and i £ [1, k — 1], we have 

h ™ — E = TTT 1 — : a ~u — T + 1-6 w? «j- (5-44) 



h-+0 h \i + 1 b — a b — a J 

Indeed, we have 

y^.(l) = / (twjYdt+ / (twj + h- — -Vj) dt+ (tWj + hvjfdt 

a rb pi 

(tw j ) i dt+ I {tw j ) i dt+ j (twjfdt 

l) J a Jb 

b + „ /•! 



+ ihw'^Vj I / f _1 - — -dt + / f _1 dt J + o(/t), as /i -»■ 0, 

/ i 5*+i _ a » ; +! 5* _ a * A 

= 2/m(1) + - — 7 — ; a ~, h 1 - &M w\~ l Vjh + o(h), as /i ->■ 0, 

\? + 1 — a 6 — a ) 



where o(/i) vanishes as h goes to zero. This proves (|5.44p . Analogously, 



x h n+l {\) - x n+1 {\) _ / k b k+l - a k+l _ b k - a k 
h-th h \ k + 1 6 — a 6 — a 



nm -n+iv^ -» + uv = 1 * » u a ^_^ + l _ 6 M |«;| fc - 2 (wJ, w ), (5.45) 
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when considering system (|1.3p . and 



lim 



x^ +1 (l) - x n+1 (l) 



h 



k b k+1 - a k+1 
k + l b-a 



b k - a k 

'-7 

b — a 






?'; 



j- 



(5.46) 



in the case of (|1.4p . Note that for all i G [1, fc] one has 

6 4 - o* 



i 6 i+1 - a i+1 



0(6 -a), 



as 6 — a — > 0, 



i + 1 b — a b — a 

for any i£{l,kj. Then, from (fOHJ) . (|5Tii|) . (f5T4"5l) . in the case ([Q]> . it follows that 

D£(w)w = U, (1 - 6 fc ) |uJ| fc - 2 (HJ, v), (1 - 6)u, (1 - b 2 )w lVl , . . . , (1 - 6 2 )«7 n v n , 

. . . , (1 - b k - l )w k - 2 v^ . . . , (1 - & fc - 1 )«5*- 2 « n > ) + 0(6 - a), 



(5.47) 



(5.48) 



as b — a —)■ 0. We next choose bo, ■ ■ ■ , bf. e]0, 1] such that bi ^ b m if i ^ m and we let v be 
any unit vector ej of the canonical basis of M. n . Then the j — th, the n + j + 1 — th . . . , the 
(k — l)n + j + 1 — th components of DE{uj)uj are 



(l,l-6 i ,(l-6fH-,...,(l-6f- 1 )^- 2 ) 



while the n + 1 — th component is (l — 6n|u;| 2 Wj- By our assumption, Wj ^ 0, and the 
following (k + 1) x (k + 1) matrix 



/ 1 1 - 6 1 - &§ 



1 " b k \ 



M(ba,h,...,b k ) 



1 1-6! 1-6? ... l-b\ 
\l 1-6* 1-6 2 ... 1-6*7 



is non singular, since 



detM(6 ,6i 



-i) fc n a* 



¥=0, 



i^m 



because of our choice of the 6j's. Thus, if we choose v = ej and each etj sufficiently close to 
bi, then (|5.48p restores k + l linearly independent vectors. In conclusion, it is possible to 
find vi, . . . , v n ,bo, . . . ,bk,ao, ■ ■ ■ , ak, such that the vectors DE{uj)uj defined by using Vj,a,i, bi 
in (|5.42p . span M +1 . This proves our claim for system (jl.3p . The proof in the case (jl.4p is 
analogous, we only need to replace (J5.45P by (|5.46p . We omit the details. □ 

We next obtain, as a corollary, a Harnack inequality which is invariant with respect to 
the Lie group G = (M fcn+2 ,o, (5a)a>o)- For every compact subset K of the unit cylinder O 
defined in (|5.40p . any positive r and any zq = (xo,yo,to) E R fcn + 2 we set 



O r (zo) = zqo S r O = {zq o 5 r ( | £ E 0}, IC r (zo) = zqo S r JC. 



(5.49) 
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We recall the definition of the end-point map E introduced in (|5.4ip , and we define the 
following sets 

I:= iuj : [0,|] ->M. n | u(a) = u with \ < \ujj\ < l,j G [l,n]|, 

I := L : [0, |] -* E n | UJj(s) = c for any s G [0, ±] U [|, ±] , 

L . (5.50) 

ujj(s) = —c for any s G [|, |] , j G [l,n] with 1 < c < v2, k 

£:={ (£?(w,0,0,0,|),i) \uel\ £:= j (£?(w, 0, 0,0, §), |) | w G /}. 

We remark that /C and K, are compact subset of the unit cylinder defined in (|5.40p . being E a 
continuous map. We also note that, if we denote by z = (x,y,t) the point (E(lo, 0, 0, 0, |), ^) 
with u) as described in I, we have xi, n = 0i, n ,x n+ i = c a n +i,fe> with: 

a»+l,* = f^8-( fe+1 ) forfEID, 

0*4-1,* = FFr8-( fc+1 ) for (gSD and k even, (5.51) 

fln+i,fc = for (|2.8p and k odd. 



Proposition 5.15 Let O C R fcn + 2 fr e i/ie unit cylinder defined in (|5.40|) . let r > and /et 
zo = (^0)2/0) *o) S M Ar+1 . If J£ is the lifted operator of \S£ in (|2.7p or (|2.8|) . iften i/iere exists 
a positive constant c such that the sets /C and /C defined in ()5.50p are compact subsets of 
Int (.2/(0,0,0) )• Moreover there exist two positive constants C^,Cp, only depending on O and 
on Jzf , suc/i iaai 

sup u < C^u^o), sup u < Cj^u(zq), 

ICr(z ) lCr(z ) 

for every positive solution u of Lu = in O t (zq). 

Proof. By the invariance of Jz? with respect to the homogeneous Lie group G, it is not 
restrictive to assume (xo,yo>^o) = (0>0, 0) an d r = 1. By Lemma 15.141 E(uS) is an interior 
point of =2/(0,0,0) f° r anv ^ G /, and for any u E I. The conclusion follows from Theorem 15.21 

□ 

6 Harnack chains and lower bounds 

In this section we build Harnack chains and we prove asymptotic lower bounds for positive 
solutions to Jzfu = 0. In the first Lemma 16.11 we capture paths that give Gaussian lower 
bounds as \xi n — £i, n | 2 > K(t — r), for suitably big K and asymptotic bounds for points x, £ 
with \x n+ i — £ n +il suitably big with respect to (t — r) 1+k ' 2 . 

Lemma [6. 21 applies when \x n+ \ — £ n +i\ is small with respect to (t — r) 1+k ' 2 . Moreover, in 
this lemma, we consider points with non degenerate components set to zero. In such case, if 
we denote (x,y,t) = (x,y,t) o (x,y,t), then 

»l,n = £l,n = 0i, n => X n+ i = X n+ i + X n+ i, (6.1) 
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that is the group law o is additive w.r.t. the (n + l)th component. In some sense this allows 
to move in the direction of the vector field [d Xl , [d Xl , • • • , [d Xl , Z] ■ ■ ■ ]] = k\d x +1 . 

V v ' 

k times 

The proof of the two lemmas is based on the lifting procedure introduced in Section 5 
and on the construction of a finite sequence of cylinders contained in the lifted domain of the 
solution. Specifically, we find points along the trajectory of the integral path introduced in 
(|5.34p . The bounds depend on the length of the Harnack chain that in turns depends on the 
slope of the trajectory. Asymptotic lower bounds are proved in Propositions 16.31 and HT 



u 



Lemma 6.1 Let ££ be the operator defined in (|2,7p or in (J2.8J) . let Ti,T,t,T 2 be such that 
Ti < r <t <T 2 andt-T < r -T 1 . Let 7 : [0, t - r] -» R n+1 x ]T 1 ,T 2 [ be a path satisfying 

n 
1 >{s) = Y j u ] Y 3 { 1 { S )) + Z{ 1 {s)) 1 j(0) = (x,t), 7 (t-r) = (e,r), (6.2) 

i=i 

for some constant vector ofM. n such that 

max I u; 7 - 1 < 2 min luA, and that (t — r) max uj,- > 4. 
36[l,n] i6[l,n] je[L,n] J 

Then there exists a positive constant C , only depending on ££ , such that: 

<■(£■, t) < exp I C((t — t) max u H + l) ) u(x, t), 
V i6[i,»] 3 J 

for every non-negative solution u to 5£u = in M™ +1 x ]Ti,T2[. 

Proof. Define the function u by setting u(x, y, t) = u(x, t) for every (x, y, t) G ]R fen+1 x ]Tx, T 2 [. 
Clearly, wis a non-negative solution to S£u = 0. Let 7 : [0, t — t] — > ]R +1 x ]Ti,T2[ be the 
solution of the Cauchy problem 

il'(s) = Ei=iWi?i(7(s)) + Zms)), se[0,t- r], 
\j{0) = {x,0,t), 

where u is the constant vector in (|6.2p . Note that, if xi n +i and xi jn +i are the first n + 1 

components of 7 and 7, respectively, then we have Xi, n +i(s) = %i,n+i(s), for every s G [0, t—r]. 

We next apply the Harnack inequalities stated in Proposition 15.151 to a suitable set of 

points z\, . . . , z m lying on 7QO, t — r]). We suppose u^ = max je j lTl ] cj|, as it is not restrictive, 

and we let m be the unique positive integer such that m — 1 < - — rf- 1 - £ 7Tl - We set s = ^^ 
and we define Zj = 7(^5) for j € [l,m]. In order to apply Proposition 15.151 we put r = v2s, 
and u = T^LO. Note that m > - — p^ 1 > 2 and, as a consequence, we have uj G J. Thus, if 
we denote by z the point (E(uj, 0,0,0, 5), |) defined in (|5.50p . we have z £ K.. Moreover, 

Zj = Zj-l O 6rZ, thus 

Zj G /Cf(zj_i), for any j G [l,m]. 
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Note that by our assumption, %<\<t-T<T-T x , hence, 0?(zj) C R kn+1 x]T 1 ,T 2 [ 



for every j E [0,m— 1]. Thus, by Proposition 15.151 there exists a constant Cp > 1 such 

that u(zj) < Cj^u(zj^i) for every j E [1,to]. In particular, being m < - — T ^ L + 1,Zq = 
(x,0,t), z m = j(t — r), we find 

u(7(t-r))<C^^ +1 u(x,0,t). 
Hence, 

u(& r) = u( 7 (t - r)) = «(7(t - r)) < CC - ^* 1 ^, 0, i) = C^^ 4 " 1 ^, t), 

and our claim follows by choosing C := log(Cp). □ 

Note that, whenever Lemma [6TT1 applies, we have x\^ n {t — r) = x\^ n + (t — t)uj / x\^ n . 
Next result gives a bound along a trajectory 7 such that xi )n (s) = for any s E [0, t — r]. 

Lemma 6.2 Let _£? be the operator defined in (j'2.7p or in (I2.8p . wii/i k even, and Ze£ a n+ i^ 
be as in (|5.5ip . Lei T\,r,t,T2 be such that T\ < r < £ < T2 and t — t < t — T\. Then there 
exists a positive constant C , only depending on Jz? , suc/i £/ia£: 

/ \ / ((t-T) 1+2/k \\ 

u ^0i )n ,x n+ i +Cn+i,Tj < exp ( C( ^ \- 1 j j n(0i, n , x n+ i , i) . 

/or every (x,t) E M n+1 x ]Ti, T 2 [,t e]T x ,t[ with £ n+1 E ]0, 2a n+ i, fc (t - r) 1+fc / 2 [, /or every 
non-negative solution u to ££u = in M. n+1 x ]T\,T2[. 

Proof. As in the proof of Lemma[67TJ we consider the function u defined as u(x, y, t) = u(x, t) 
for every (x,y,t) E M fcn+1 x ]Ti,T2[. Let m be the unique positive integer such that 

/ \ 2/k 

ro -l<(2(t-r)) 1+a/ * ^±M < m . ( 6 .3) 

Let r = \ and let c = ( ) m 1/2+1/fc • Note that from our assumption 

V m \a n+1 , k j ( 2 (i_ T )) 

< Cn+i < 2a n+ i j fc(t — T) 1+fc ' 2 it follows that m > 2, hence 1 < c < \/2, Then, if we denote 
by z the point (E{uj, 0, 0, 0, |), |) defined in (|5.50p . with c as above, we have z E /C. 

Let zo = (0i jn ,x n+ i,0, £), and let Zj = Zj-i o 5r~z, for j E [l,m]. By our assumption we 
also have r 2 < r - Ti, then Oyr(^j) C M fcn+1 x]Ti,T2[, for every j E [0,m]. Thus Proposition 
15.151 yields u(z m ) < C^u(zo). According with (|6.ip . and with our choice of m,c, and r, the 



first n + 1 components of z m are ( 0i iTl , x n+ \ + £ n+ i j . Then 

u(0i >n ,x n+ i +x n+ i,r) = u(z m ) < C~u(z ) = C~ u(0i, n ,x n +i,t), 
and the conclusion follows from (16.31). □ 
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Proposition 6.3 Let Jzf be the operator defined in (|2.7p and let k be a positive even integer. 
Let u : R n+1 x]Ti,T2[— > R 6e a non-negative solution to ££u = 0, and let t, r G R 6e suc/i 
t/iai Ti < r < t < Ti, and t — r < 2(r — Ti). TTien i/iere exists a positive constant C\, only 
depending on J£ , such that 

i) foranyx,££R n+l such that £ n+i - x n+1 > 15 fc ((i-r)(|xi in | fc + |a,n| fc ) + n fc / 2 (t-r) 1+fc / 2 ) 
we have 



/ _ \ 2 / fe 
ufc, r) < exp \C X ^ — + _ r)1+2/fc + 1 jj 



2/fc 

u(x,i); 



nj /or any z,£ G M n+1 suc/i i/iai < £ n+1 - x n+i < ^TT) Xl >«l fe + l£i,n| fc ) + 8 k+i {k+1) {t ~ 
r) 1+fc / 2 we have 

u(f,r) <exp Ci + — — ^ + 1) )u(x,t). 

Proof of (i). We divide the proof into two steps. In the first one we find a path 7 : [0, ^-] — > 
j^fcn+2 ^at steers xi >n to £i )n . In the second step another path 71 + 72 steers the (re + l)th 
component x n+1 (^p) of 7 (^) to £ n+i . 

Step 1: If w = tz^{xi jTI — £i, n ) satisfies the assumptions of Lemma |6.1| then one read- 
ily gets zi, n (^r) = Ci,n and u ( 7 (£=!)) < exp(C[2 |a:i -" t ~^" 12 + l])«(ar,t) and the first 
step is achieved. If this is not the case, we rely on the following construction. Set K = 
max je j ln j /,_ and M := max{3-ftT, |}. For every j G [1, re], we set 

AM . A Cj- Xj AM 



, ujj := 4 ;= 



so that 



8 M 16 M i-r„ 9 £-t_ 9 

Consider now the path 7 associated to (|6.2p with w(s) = wl se r t^i +wI se rt ZI T t-ri for which 
x i,n {^y~) = £l,n- The assumptions of Lemma 16. II are clearly satisfied. Hence: 



(7 {*?)) < exp (c> C^; _^' w| + l)) <x,t), 

j( t ^ L ) = Ul,n,^+1+ / Xl,n+ / u(u)du 



1 t + A (6 " 4) 



for some positive constant C . By a plain change of variable in the above integral we find 

t—T t — T 

r 4 , r 4 , 

^n+i (^jr) - x n+ i = I \x\ jU + su\ ds+ \£ 1>n - sQ\ ds. (6.5) 

Jo Jo 



'0 
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Note that, for s £ [0, ^-r 1 ] and j 6 [1, n], we have 



\xj + suij\ < \xj\ + -4=M < \xj\ + 3 I \xj - 01 + ^V") < 4|xj| + 3|0l + 3^ 



t-T 



\ x j-ij\ 



10 - *%| < 101 + 7/fe M + ^t=T ^ 4|x,| + 5|0I +3 



2 ' 



(6.6) 



thus, from the elementary inequality (a 2 + b 2 + c 2 ) fc / 2 < 3 fc / 2 1 {a k + 6 fc + c fc ), we get 

( t - T \ ^. ^ t-rok-l (r ) 2k+l\ rr , \k j_ (ok j_rh\\(: \k , n /^n*-/^ 

^n+1 (,— J - ^n+1 < — 3 12 |»i jn | +{6 +b J|£l,„| + 2 1 2 fc M 

Recalling that £ n+i - x n+ i > 15 fc ((i - r)(|xi in | fc + |0,n| fc ) + n fc/2 (i - r) 1+fc / 2 ), we find 



t-T 



£n+i - x n+1 (t=f) > 15 k (\ Xl , n \ k + |ei,„| fc + n k / 2 (t - rf' 2 ). 



(6.7) 



We next prove a similar lower bound for a; n +i (^J-j — #n+i- To this aim, we note that, if 
\xi,n\ > ^-jp|w|, then |xi )n + suj\ > sjNi.nl f° r every s S [0, £j#] . Analogously, if |xi )R | < 
^p|u;|, then \xi lTl + su;| > ^jjp|55| for every s£ [^(* ~~ r ); ^r~] • The same remark holds if we 
replace xi <n and ui by 0,n and u), respectively. As a consequence we find, 



/*__\ t — T ( \Xl 

x n+ i - x n+1 {-*-) > —r^r- I max 



16 



2 k 



so that, in particular, 



t-T, 
~~16~ 



t-T 



\U\ 



+ max 



ic 



l,n 



2 A; 



t-T 

16 ' 



X n+ i (^) - X n+ i > -r-r^ (Jxi,„| + |0,n 



(6. 



Step £: We next denote by u the vector in W 1 such that cjj = 1 if > 0, Wj = — 1 if 
< 0, for j £ [l,ra], and we fix a real parameter 6 > , 4 _ , that will be specified later. We 
consider the path 71 : [^, |(£ — r)] — ► M n+1 x]Ti,r 2 [, starting from 7 (^ 2: ) and defined as 
in (pT2|) with w = 6uJ, then the path 72 : [|(t-r),i-r] -> W n+1 x]T 1 ,T 2 [, starting from 
71 (j(t — r)) and defined by setting uj = —bio. Prom Lemma 16. II it then follows 



u(£l,n,<P(b),T) < 



exp( 2 c(^^ + l))«( 7 (*?)), 



(6.9) 



where 



<p(b) = x n+1 (±=f) + 2 / \£ lin + sbui\ K ds 



t-T 

4 



is an increasing continuous function of b 6 
that 



/t-T 



+ OO 



An elementary computation shows 



^■- i^.^o**-^., .* , & fe n fe / 2 (t - r)^ 1 



|6,n + «M d« < 2*— -|6,nr + 



k + 1 2 fc + 2 
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then 



V I -7^=) < s.H-1 (V) + 2^|6,n| fc + ^^(t - r) 1+fe / 2 < £ n+1 , 



by ()6.7p . On the other hand we have 

b k u _ T \k+i n k/2 
?(&) > z„+i (*?) + ^-j-^ ^fc+i ^+°°' as 6 ^+00. (6.10) 

Hence, there exists a unique b > j== such that (p(b) = £ n +l- Moreover from (|6.10p it also 
follows that 

b < 2(fc + l) 1/fc (&) W/k (Cn+i - x n+1 (^)) 1/fc n" 1 /*, 

which, together with (|6.8p , gives 

1 4-1 Ih / + \ Iffe 

b<(k + if' k (£) ^ n+1 - x„, +1 - ^J (|x ljn | fc + iei,„i fe ) j . 

Eventually, equation (|6.9|) yields 

n(e,r) < exp (^ (^ (t - r )i+2/fc + l ) ) u fr (V)) ■ 

with Ci = 4 1+1 ' fe (/c + l) 2 ' k C '. The above inequality, with (|6.4[) . proves the claim (i). 

Proof of (ii). We prove our claim by applying Lemma f6.2l in a suitable interval [r+foj t— 1±] C 
[r, t], and Lemma 16. II in the remaining intervals [t — t\,t] and [t,t + £ 2 ]- We first suppose 
that xi >n / 0, £ lin 7^ 0, we set 

. J max j£[l,nJ x j Cn+1 - ^n+1 * ~ T \ 
t\ = mm ' v 



I max 
ti = min < - 



\xi,n\ k 



■ie[l,»] £j Cn+l - £n+l * - T 



16 '* 



3l,n 

and we consider the paths 

, n /7, s \ (s - £i) A:+1 + *i +1 , lfe \ rn , 

71 (s) = I ^1 - y) Xi ^ Xn+l H rfc + ii^ — ' ' * ~ s 1 ' s G 

fe+i + fc+i 



s . s K+1 - t 



72(s) = -gl,n,gn+l + ,, , ,A |$l,nr,r + t 2 -8 , «£ [0,t 2 ]- 



£ 2 -™^ ' (fc + l)f 



We next proceed assuming that max je | ln j \xj\ < 2minj e | ln j \xi\ and max je | ln j |£j| < 
2 mirijgnjj] |£j|. In this case we apply Lemma IBTTl in the interval [0,£i] with uj = —j- x l,ni so 
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that we have t\ maxj e [ lin ] u, > 4 and max je j ln ] \ujj\ < 2 min !e [ lin ] |wj|. We find 



u(£, t) < exp I C I max 






Sn+1 ^n+1 

with 



' 3 T i 4' 4 } + 1 )) u(72(0)): 



7i(*i)= (0,x„ + i + ] # T |xi jn |V-tiJ, 72(0) = (^0,e„+i-Ff T |6 i nr,r + t2 

If maxj g | 1)n j |xj| > 2min ie j lri ] |xj|, we rely on the argument used at the beginning of the proof 
of (%). Specifically, we set M := max{3maxj e | ln j ^U., |}, and uij := 2-4=, tDj := pXj + 2—== 
for every j £ [l,n], then we consider the path 7 associated to (|6.2p with u)(s) = £31 rn *ii — 
SI r ti . ,. Also in this case we get (|6.11[) . with some bigger constant C. Aiming to simplify 

our exposition we omit the details of the proof. 

We next conclude the proof by using Lemma llOl We set £ n+ i = ^ n+ i — x n+ i — -^^\x\ in \ k — 

Efr|fl, n | fc , and we recall that £ n+ i - x n+1 < ^pjy (|xi,„|* + \£i, n \ k ) + 8 ^(k+i) (* ~ r) 1+k / 2 . 
Thus 

L_l < ( t _ tl ) _ ( r + t2 ) < t _ T) _1 ( ?n+1 _ Xn+1 ) < ^ n+1 < _^__( t _ r )i+*/2. 

(6.12) 
Then, from Lemma 16.21 we get 

« (75,(0)) < exp (c ( (t "' 1 ( ~' 2 i) l / I )1+2/fc + X )) "(^0) 
From inequalities (|6.11|) and (|6.13p it follows that 

/ /| T1 |*+2 i I* |fc+2 /y_ \l+2/fc I | 2 + |A I 2 

u(f , r < exp C + — — - + + 1 u(x, t 

V V Wl-^n+1 {in+l-X n+ iYl k t-T 



(6.13) 



7 v ) •> 



for some positive constant C only depending on C and on k. Note that the last term in the 
above expression is bounded by the first one. Indeed, the inequality 

| Tl | 2 -l|A I 2 9 In | fc + 2 + |A l fc + 2 £• 

1^1, n\ r |?l,n| ^ ^ l x l,n| t lsl,ra| , "> 



t-r ~ k + 2 (i_ r )i+fe/2 ' fc + 2 

combined with (|6.12p . gives 

|X1,„| 2 + |ei,n| 2 < 2 |si,n| fc+2 + |4~l,n|* +2 + fc 



t _ r - 4*+ 1 (fc + 2)(fc-l) Cn+i-x n+1 k + 2' 
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This concludes the proof of (ii) when xi >n 7^ 0, and £i )7l 7^ 0. 

If x\ t n = 0, we simply omit the construction of 71, and we rely on 72 and on the application 
of Lemma 16, 2 1 in the interval [r + t 2 , t]. Analogously, if £i )n = 0, we avoid the construction of 
72. This concludes the proof. □ 



Proposition 6.4 Let «Sf be the operator defined in (|2.8p and let k be a positive integer. Let 
u : M. n+1 x]Ti,T2[— > R k o non-negative solution to ££u = 0, and let t, r € R be such that 
T\ < t < t < T2, and t — T < 2(r — Ti). TTien £/iere exists a positive constant C\, only 
depending on ££ , such that 



i) if k is even, then for any x, £ £ M ra+1 suc/i i/iai £n+i — £n+i > 15 (t — r) X/?=i pHnK 
n (|) (t - r) 1+fc/2 we We 



u(£, r) < exp ( d [ — \_ T 1- 



(S n+1 - x n+l - ^ E -=i(4 + #)(« - r)?l k 

+ (t- T y+^ — + ! ) " ( "''- /) - 



mJ ?/fc is even, then for any x,£ £ R n+ such that < £ n+ i—x n+ i < 2 tk+i) 2j=l l 2 ^ + £f ) + 



(t _ r) l+fc/2 



,,., |xi, B | fc+2 + |gl,n| fc+2 , (t-T) 1 ^ , ,, , , 

u(£,t) <exp Ci + — r^yr + 1 u(M); 



Hi) if k is odd, we have 

u(tr)<e^(c 1 Xl ' n ~_^ nl \ 

+ ( t - T jl + 2/ fc + 1JJ«(^*)- 

Proof. The proof of (%) and (ii) is analogous to that of Proposition 16.31 The unique modifi- 
cation is due to the fact that here we have 



~*ziZ n „LzZ n 

Wi (^r) - x n+ i = Y] (xj + su>j) k ds+ Y] (£j - su)j) k ds. (6.14) 

Jo ~1 Jo ~T 



instead of (J6.5J) . From (|6.6|) , by the same argument used in the proof of (J6.7J) and (j6.8|) , we 
obtain 

t — T \-~\ ( u .j,\ , /•+_— \ t — T 



£ (4 + tf) < ^+1 (*?) -*„ + i < ^15 fc £ (l^l* + |fc|*) +2n (§)* (t-r)^ 1 . 



2*; +4 
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We omit the other details of the proof of (i) and (ii). 

The proof of (Hi) follows from the same argument used in the proof of (i). Note that, as 
k is odd, the function 



« - — - it, 

<p(b) = x n+1 (^) + 2 / ' y^ixj + bWjfds 

JO .-_! 



defined for any b E R is surjective, then in this case Lemma 16. II is sufficient to conclude the 
proof. □ 

Final derivation of the estimates 

Proof of Theorem \2.1\ It follows from the bounds proved in Sections 14.51 and Consider 
first equation (|1.3|) . 

The upper bound of (i) is given in (|4.38[) . In order to prove the lower bound, we fix a 
constant vector to £ M. n ,a G]0, 1[, and we set 

r f a2t k 

Xl,n=€l,n-CtVtU), X n+ l=Cn+l- £l,n ~ ^W ds. (6.15) 

According with Remark 15.131 we have 

p(a 2 t, x, = T(x, a\ g, 0) = - -^r , C = r(w,-jgj,l) (6.16) 

(a 2 t) + 2 v ' 

For our purpose, we choose here u = (1, . . . , 1) € R n , and a = l/v2. Note that the constant 
C is strictly positive, being uj ^ 0. Also note that 



< £n +1 - x n+1 < 2 fc a 2 (t|ei,„| fc + a fc t 1+fc />| fe ) . (6.17) 

We then apply Proposition 16.31 (i), and we find 

/, /-\ \ I ri I \ x l,n ~ x l,n\ 

p{t, x, £) > exp I - d I h 



2/fc 



I 2-n+l ^n+1 ofc+5 ( P'l.n ~r P^l.n J 

V '—^ + l))p(t/2,x,0- 



t l+2/k 

The lower bound (i) thus follows from the inequalities \x\, n — £l,n| < vt/2 M and £, n +i > x-n+i- 
The upper bound of (ii) is equation (|4.37p . the lower bound is given in Lemma I4T71 
The upper bound of (Hi) is given in (|4.39p . In order to prove the lower bound we rely 

on Proposition 16.31 (ii). In order to satisfy its hypothesis, we choose a G]0, 1/V2] such that 

in+i ~ x n +i < \{in+i - x n+ i). Note that 



< U+i ~ xn+i < Ta l [t\x hn f + a*t 1+ */>|* ) , (6.18) 
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then it is sufficient to choose a such that 






1 . Cx,k I \ x l,n\ , 

—7^- > max 7— , 1 

a 2 t ~ t I - ■-*-' 



Using our assumption |£ n +l — ^n+l| < -ftT£ 1+fc ' 2 we find 

\k 

for some positive constant Cx fe depending on K and fe. Then (I6.16P gives 

_ (fc + 2)(n+fc + 2) 

rV ■ / l^r, l fc + 2 \ 2fe 

p(a t, X, > — p^r 1 + y 

t — 2 — V l?n+l - Xn+1 

With this choice of a Proposition 16,31 (ii) then gives 

_ (fc + 2)(n + fc + 2) 

Ck k ( 15? 1 \ k+2 ^ 2fc 

p(t,X,0 >^TSTF 1 + 17 != — 



^+2 i I™, |fc+2 

exp - Ci '— '■ h -pz r^-pr + 1 



|xi,„|— + |xi, n |— , ((l-a 2 )t) 1+2 / fc 



X n+ i - X n+ i {Xn+1 ~ X n +l) 2/k 

and our lower bound follows from the inequalities t/2 < (1 — a 2 )t < t, l/2(£ n+ i — x n +i) < 

Xn+1 ~ Xn+1 < £n+l ~ X n +1 and |xi, n - £l, n | < V^/2 |w|. 

When considering equation (|1.4p . the lower bounds for points (%) and (mj directly follow 
from 16.41 The proof of the remaining bounds can be done by the same arguments used for 
equation (jl.3p . □ 

Remark 6.5 We can assume by symmetry that w.l.o.g. |£i n | > |xi,n|- In this case, observe 
from equation (|6.17|) that i/|£i n | > Kt 1 ' 2 for K large enough, then we can derive from Lemma 
\4-7\ that the Gaussian diagonal regime holds for the lower bound. From the proof leading to 
(|4.38p . this means that the non- exponential estimates in case i) could be alternatively rewritten 
changing the i - "- ~ + > term into 



r n/2 t -3/2 ({ | Xijn | fe -l + | 6n |fe-l} V t k ~^)-\ (6.19) 

that emphasizes the regime transition depending on the magnitude of the non- degenerate com- 
ponents w.r.t. to their characteristic time- scale. From the above computations, the expression 
in (|6.19p can also substitute the non- exponential term in case Hi). 
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